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Ein Verfahren 
zur Stabilitatsfrage bei Matrizen-Eigenwertproblemen 


Von Hans-RuDOLF ScHwarz, Ziirich?*) 


"9 1, Einleitung und Problemstellung 


Zahlreiche technische Fragen fiihren auf das mathematische Problem, fiir. 
das Matrizen-Eigenwertproblem (A —AE) x = 0 die Anzahl der Eigenwerte - 
“mit positivem Realteil zu bestimmen, wobei die Eigenwerte selbst meistens gar 
nicht interessieren. Diese Problemstellung nennt man kurz die Stabilitatsfrage. 
_Um diese zu beantworten, sind verschiedene klassische. Kriterien bekannt, - 
~ welche aber alle vom charakteristischen Polynom P(A) ausgehen, also bedingen, 
_ dass dasselbe explizit ermittelt wird. So erwahne ich die Kriterien von Routs{[1]2) 
und Hurwitz [2], welche in ihrer urspriinglichen Formulierung auf Polynome ‘ 
mit reellen Koeffizienten Anwendung finden, im iibrigen aber gleichbedeutend . 
sind. Das Ortskurvenkriterium von Nyquist [3] bestimmt auf graphischem 
- Weg die Anzahl der Nullstellen der charakteristischen Gleichung P(A) = 0 mit 
‘negativem Realteil. Eine unmittelbare Folge davon ist das Lagen- oder Liik- 
_kenkriterium [4], welches aussagt, dass dann und nur dann samtliche Null- 
stellen von P(A) = 0 negative Realteile aufweisen, wenn die Wurzeln von Real- 
und Imaginarteil von P(t y) = Re(y) + 7 Im(y) reell sind und sich gleichzeitig 
_ gegenseitig trennen. Das bekannteste numerische Verfahren zur Beantwortung 
der Stabilitaétsfrage diirfte wohl die Methode der Sturmschen Kette [5] sein, 
_wonach aus dem Real- und Imaginarteil von P(¢ y) = Re(y) +7 Im(y) durch 
_fortgesetzte Division mit Rest eine Folge von Polynomen absteigenden Grades 
-_-gebildet wird. Die Anzahl der Eigenwerte in der linken komplexen Halbebene 
Rez < 0 lasst sich dann aus dem Unterschied der Zeichenwechsel in der Poly- 
; nomfolge fiir y= —oo und y= +00 angeben. Damit dquivalent ist das. 
-Stabilitatskriterium von WALL [6], welches die Aussage in die Form eines 
. Kettenbruches kleidet. Daneben finden sich in der Literatur einerseits noch 
_verschiedene numerische Verfahren, welche vermittels eines Reduktionspro- 
_zesses aus dem charakteristischen Polynom eine Folge von Polynomen bestim-: 
: men, aus denen die Frage beantwortet werden kann, und anderseits auch 
| einige graphische Methoden. 
J *) Institut fiir angewandte Mathematik der ETH. 
: 2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 500. 
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Auf der andern Seite fiihrt jedes gangbare numerische Verfahren zur expli- 
ziten Entwicklung des charakteristischen Polynoms von (A = AE) «x=0 auf 
eine Folge von Polynomen, deren Grad sich bei jedem Schritt um 1 erhoht 
und die mit dem charakteristischen Polynom endigt. Ich verweise auf das Ver- 
fahren von WEBER-VOETTER [7], auf die Methode von HESSENBERG [8] und den 
Biorthogonalisierungsprozess von LANczos [9], welche alle tiber eine Reihe von 
Polynomen aufsteigenden Grades zum charakteristischen Polynom fithren. 

Um dann die allein interessierende Stabilitatsfrage beantworten zu kénnen, 
muss eine Polynomfolge absteigenden Grades ermittelt werden. Um den Weg 
iiber eine aufsteigende und dann iiber eine absteigende Kette von Polynomen 
zu vermeiden, in welcher das charakteristische Polynom als solches von héch- 
stem Grad meist gar nicht interessiert, gab mir Herr Professor Dr. E. STIEFEL” 
(ETH) die Anregung, folgende Frage zu untersuchen: 

Problem: Man berechne das charakteristische Polynom P(A) aus der Matrix A 
tiber eine solche aufsteigende Kette von Polynomen, aus welcher schon selbst die 
Stabilitatsfrage beantwortet werden kann. 

Dieser Wunsch fiihrte zu einem Verfahren, unter Umgehung des charak-— 
teristischen Polynoms die gegebene Matrix A mit Hilfe von elementaren Trans- 
formationen auf eine Normalform zu bringen, deren Elemente allein schon das 
Stabilitatsproblem lésen. Die beschriebene Methode weist den Vorteil auf, dass 
die anzuwendenden Rechenoperationen von grosser Einheitlichkeit sind, was 
das Programmieren auf Rechenautomaten erleichtert. Bei der klassischen 
Methode der Sturmschen Kette stiitzt sich die ganze Entscheidung der Stabili- 
tat auf eine langere Rechnung, welche vom charakteristischen Polynom aus- 
geht. Die Koeffizienten dieses Polynoms enthalten sicher weniger Information 
als die urspriinglich gegebenen Matrixelemente. Es ist zu hoffen, dass unsere 
Methode, die die Matrixelemente in mehr direkter Weise verwendet, in manchen 
Fallen numerisch stabiler verlauft. Ob diese Hoffnung berechtigt ist, kénnen 
allerdings erst umfangreichere Versuche entscheiden. 

An dieser Stelle méchte ich den Herren Professoren Dr. E. STIEFEL und 
Dr. H. RUTISHAUSER, die mit ihren wertvollen Hinweisen zum Gelingen dieser 
Arbeit beigetragen haben, meinen besten Dank aussprechen. 


} 

: 

I. Das Stabilitétskriterium : 

; 

2. Das Kettenbruchkriterium von WALL [6] ; 


Es sei P(e) SF Og Se Cy ig ee ee (12+, ein Polynom n-ten 
Grades mit komplexen Koeffizienten c, = p, +79, Aus diesem bilde man 
Q(z) = [P(z) — (—1)" P(—2)]/2, worin P(z) das Polynom mit den konjugier 
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komplexen Koeffizienten ¢, = p, — i qx bedeutet. Dann ist 
Q(z) = a gn-l a t Cate gn-2 = Powe gn-3 are 1 Pyne gn-4 Se 


und der Quotient Q(z)/P(z) hat im allgemeinen eine Kettenbruchentwicklung 
von folgender Form mit in z linearen Teilnennern: 


Q( 


P( 


2) 1 | 1 ] 1 | 1 | 
Dole +s, 1 is aes ot l7s2 + Sy Spee ia ad ey 

Satz 1 (Kriterium von Watt): Existiert die Kettenbruchentwicklung (2.1) 
von Q(z)/P(z), so sind darin die Grissen 7, (k = 1, 2,...,n) reell und von Null 
verschieden, die s; (k = 1, 2,..., ) rein imagindr oder null. Sind dann unter den 
%, mM positiv und (n — m) negativ, so haben m der Wurzeln von P(z) = 0 negative 
und (n — m) positive Realteile. 


Die Teilnenner in (2.1) ergeben sich durch fortgesetzte Division mit Rest 
nach dem euklidischen Algorithmus fiir P(z) und Q(z). So ist speziell fiir den 
ersten Teilnenner 


P(z) = (7, 2 + 8, + 1) Q(z) + Ril) (2.2) 
und nach Subtraktion von Q(z) auf beiden Seiten 


P(z) — Q(2) = (41 2 + 81) QO) + Ril) - 24) 5 


Da der Rest R,(z) in (2.2) und (2.3) gleich ist, lautet die Kettenbruchentwick- 
lung von Q(z)/[P(z) — Q(2)] 


O(z) iy L 1 | 1 PEGE 1 | 


P(z) — Q(z) Lege 385 |%> 2+ So te es We Soras 


wy 


(2.4) ; 


mit denselben 7;,, s;, (k = 1, 2,...,) wie in (2.1). Demnach gilt der zu Satz 1 
aquivalente 


Satz 2: Existiert die Kettenbruchentwicklung (2.4) fiir den Quotienten 
Q(z)/[P(z2) — Q(2)], so sind darin die Grossen 7, (k = 1, 2,...,) reell und von 
Null verschieden, die s;, (k = 1, 2,..., 1) rein imagindr oder null. Die Anzahl der 
positiven 1, ist dann gleich der Anzahl der Wurzeln von P(z) = 0 mit negativem 


Realteil. +e 
Die Teilnenner (7; 2+ s,) in (2.4) berechnen sich nach dem euklidischen 


Algorithmus fiir f(z) = P(z) — Q(z) und f,(2) = Q(z). Es ist 
folz) = P(e) — Ql) 
= gn 1 Gay gn-l 5 Pin. 6 gn-2 al. 1 Cues gn-3 aie Pn—a gnu-4 Jo abo, (2.5) 


(2) = Ole) = Boa 24 i dnt + Pag OH Gdn eben, 


A 
i ee 
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und es sei allgemein 
n—k 
: = y _yn—k-j 
fe (2) = ay, k ie an Vy, b+4 gO lage aiige Ap, n—1 * is ann = Bn, +4 z 
, f =v 


CP = OTL paola 
Dann bildet man 


fo(2) = (11 2 + $1) f(z) + fol) » | 
f(z) = (722 + So) fol2) + f5(2) » 


site| tertet (s) o("61 OI 6ice/ 581) Bi 01/9 V0.8), OL eR = 0.80. ee ee ae kee ae 


Se. al 6) © 16).6y,0; a 0) 018) 6. (Bite) ip (e| 1e) @ | @) 10] (0 (e.e 16) 6) ane a4 


Die Gréssen 7, Sz, Gi, nag (R= 1, 2,...,237= 0, 1,-..,0— FR) lassen sich tekum- 
siv gemiss den nachstehenden Formeln berechnen: 


Ay_4, k-1 + Ur, = Vk > (Ap-1,% — "re @n, ei) * Ann = Sx (k= 1, 2).2598) 
ei, k+1+5 — Ske %,n+14+5 — "he Ve, e+24+5 = C41, R+I4G ‘ (2.7) 
(Red, 2h 154 = 0,1, nt — a a ee 


Die Anfangsbedingungen dazu lauten nach (2.5): 


Go = 1, Gy =1Gn1, %2—=Pn-2, % = In-8; Goa = Pn—a,-- 


43 =Pnas> %2=t9n-2, %N3=Pn3, UWs=tIn4,---- 


3. Herleitung des Kriteriums von WALL aus dem Verfahren der 
Sturmschen Kette 


Betrachten wir die Methode der Sturmschen Kette [5] zur Beantwortung 
der Stabilitatsfrage: Danach setzt man in P(z) z = 7 y ein, trennt in Real- und 
Imaginarteil P(t y) = Re(y) + i Im(y) und beginnt fir gerades n mit 


Fy(y) = Rely), Fy(y) = —Im(y), 


bzw. fiir ungerades m mit F(y) = Im(y), A(y) = Re(y). Dann lasst sich im 
allgemeinen eine Folge von Polynomen F,,(y) vom wirklichen Grad » — k und 
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mit reellen Koeffizienten A, ,,; mit A, + 0 wie folgt bilden: 


F,(y) = Ax, x nF Are, ng VRE $f wee H yy Wire Aga 


n—k 


= om Ay 45 ‘aa e= 4 (k = (0) i wee y n) ; 


Fo(y) = (Ri y + S,) Ky) — K(y) , 
Fi(y) = (Rp y+ Ss) Biv) — K(y) , (3.1) 


Cree Was &% eee a's tye «6% 
STG G NO eager e S18 10 Whe ap ele eke 
al e6h eS <0 lara) BP 6 ace als 


F,-4(y) a (R,, er Sn) F,(y) - 
Fiir das folgende treffen wir die 


Voraussetzung: Die Sturmsche Kette fiir das Polynom P(z) existiert in der 
Form (3.1), wobet die Quotienten (R;,y + S,) bet der fortgesetzten Division mit 
Rest linear seien. 


Dann gilt bekanntlich 


Satz 3: Evxistiert die Kette (3.1), so ist die Anzahl N der Nullstellen von 
P(z) = 0 mit negativem Realteil N = (n +1 —m)/2, wober (lL — m) der Unter- 
schied in der Zahl der Zeichenwechsel in der Sturmschen Kette Fy(y), Fy(y),.--; 
F,_1(y), aly) fiir y = —0co und y = +00 1st. 


Darin ist / die Gesamtzahl der Verluste und m die Gesamtzahl der Gewinne 

an Zeichenwechseln in der Kette (3.1) beim Durchlaufen der imaginaren Achse 

von —ioco nach +ioo. Diese Werte lassen sich aus den Koeffizienten A,, , 

bzw. aus den Quotienten Ry = Ay_y,4-1/Ax,, ohne weiteres angeben. Betrachtet 

man F,_,(y) und F,(y), so hat man beztiglich dieser beiden Glieder in der Kette 
_hinsichtlich Zeichenwechsel folgende Situation: 


i fiir 'y = —co 1 Zeichenwechsel 


| 1 Veri 
fiir y = +co 0 Zeichenwechsel { — 


sgn R;, = 
oe ea 1 J fiir y = —oo 0 Zeichenwechsel 


= 1 Gewinn . 
| fiir y= +00 1 Zeichenwechsel | ae 


-Somit ist die Gesamtzahl der Verluste 1 an Zeichenwechseln gleich der Anzahl 
der positiven R;,, und die Zahl der Gewinne m gleich der Anzahl der negativen 
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R,. Wenn F,(y) und F,(y) teilerfremd sind, die fortgesetzte Division mit Rest — 
wirklich mit F,(y) = const + 0 endigt, so ist m = 1 + m und demnach 


alae eis 2, 


Nise 2 


also folgt aus Satz 3: 


Satz 4: Die Anzahl N der Nullstellen von P(z) = 0 mit negativem Realteil ist 
gleich der Anzahl der positiven Rj, in der Kette (3.1). 


Fiir die Polynome Fj(y) und /(y) ergeben sich: 
1. Fall: n= 4k oder n= 4k+ 1 (k ganz) 


Foy) = 9" + na VO = Bag I = (nag Pe ee 
Fiy)=! Pua y+ Gna V* = Png Wed ee 
2. Fall:n=4k+ 2 0dern=4k + 3 (k ganz) 
F#(9) <9" — Gna 9"? + Pua 2+ Ing" — Pray — $$ = = Rly), 
FR)= pay Gy te 


Da F*(y) = —K(y) und F*(y) = —A(y) ist, so unterscheiden sich die 
Ketten F*(y) und F,(y) nur durch das Vorzeichen; es ist F*(y) = —F,(y), 
aber R* = R, (k= 1, 2,..., 7). Deshalb kann die Zerlegung von P(z) in Fj(y) 
und /,(y) unabhangig von in der ersten Form (3.2) geschehen. Fir die Koef- 
fizienten R;,, S,, Ag,ny5 (Rk = 1, 2,...,”; 7=0,1,...,—k) bestehen die fol- 7 
genden Rekursionsformeln: 


Agora Age = Ry, (ea, Ax, 241) > Ann = Sx 


(eerie, eos) t 
(3.3) 


== ” a, 4 
Aya, eei¢s +55 Ap eet + Re An esors = Ags, ray 


(k= 1,2... 0 = 13 7 2051) 2S ait Ay nga oO) 
Die Anfangsbedingungen dazu lauten nach wae 


Ay =1, Ay = Gn» Aggy = —Py 9; “Ay = —n-g Aes = Prous 
(3.4) 
Ayn =fo1 Ag= Qn» Aig =—Py-3, Aw= —~|n—4) +++ 


Nach diesen Vorbereitungen kann nun gezeigt werden, dass Satz 2 und 
damit Satz 1 eine Folge von Satz 4 ist. | 


Voraussetzung: Die Sturmsche Kette in der Form (3.1) fiir P(z) existiert. 
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Behauptung: Die Kettenbruchentwicklun (2.4) v 
stiert, und es gilt § (4.4) von Q(z) /[P(z) — Q(z)] exi- 


Gy n+9 =U Apnss (R=015..., 0; T= O, leet oe By 


Dy 
r, = Ry, S. = oS, (R= 1,2; at, m). z 


Bewets 


iY nduktionsannahme : Die Behauptung (3.5) sei richtig fiir k <m =1. 
Induktionsschluss: Gemass Rekursionsformel (2.7) ist fir k = m: 


m+, m4i1+5 = %m—1, m4143 — Sm Om, m4i45 — 1m Qn,mvers (OSG Sn-—m— ol) 
Nach Induktionsvoraussetzung ist dies einerseits 


cree ee 4d Ay 
U ae m+14+79 U on ait ‘a m+1+5 faiee gira An, M+245 > 


ry a ca. m+1+5 “_- ee ae 7 Kan AD meses) ) 
und anderseits infolge (3.3) 
Amis, mitts = Anis mira (f= 0,1,...,%—m—1). 
Weiter schliesst man analog unter Verwendung des letzten Resultates fiirj = 0,1: 


Tmo = Imm: m4, m1 = r. tee : An+1, m1 = Poni , 
Spa = (Gm, mer — "ng. msi, m+2) ¢ m1, m41 » 


= (1 re mL Raw U Agia, ma) > Anes, maa = 4 Sse 


Induktionsverankerung: Fir k = 0 und k= 1 liest man die Beziehungen 
Ay; = 1 Ags, 44; =U Ay,,; direkt aus den Anfangsbedingungen (2.8) und 
(3.4) ab, und es folgt daraus, dass fiir k = 1 


y= Mp 2 yy = Ay: An = R,, 
$1 = (4g, — 71 %2) ? 41 = (@ Agr — Ry t Aj) : Ay = 7 Sy 


ist. Nach Voraussetzung iiber die Existenz der Sturmschen Kette (3.1) sind 
die R, + 0(k = 1, 2,..., ”), und deshalb sind auch die ry = 0\(2 = 1, 2,565.98): 
das heisst, die Kettenbruchentwicklungen (2.4) und (2.1) existieren und die 7, 
sind reell und von Null verschieden. Da die S, (k = 1, 2,..., ) reell sind aber 
eventuell verschwinden kénnen, sind die s,= 7S, (k = 1, 2,...,) rein ima- 
ginar oder auch null. Infolge der Gleichheit 7, = Ry decken sich die Aussagen 


.) 
x | 


t 
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iiber die Anzahl der Wurzeln von P(z) = 0 mit negativem Realteil in den 
Satzen 1 und 4, so dass damit der Beweis von Satz 1 mit Hilfe des Satzes iiber 
die Sturmsche Kette erbracht ist. 


4, Folgerung aus dem Kettenbruchkriterium 


V oraussetzung: Die Entwicklung von Q(z)/P(z) in einen Kettenbruch (2.1) 
existiert. : 
Folgerungen: Nach Satz 1 sind die Gréssen 7, (k = 1, 2,...,) von Null 
verschieden. Deshalb lasst sich auf den Kettenbruch (2.1) eine Aquivalenz- 
transformation ausfiihren, indem der k-te Teilbruch mit —1/r, erweitert wird: 


aN gata 


Q(z) —1/r, | 
see fp one 4.1 
Pl) VEAP SoH : aah Tae (4.1) 
Va V4 Ts 1s ve 
Setzt man dann 
=i 
= ity : = ay amo foe 
Vy Vyp—-1 an (4 2) 
= Bab, (b=4, 2 ee | 
Vn, 


so sind die a;,(k=1,2,...,m) reell und von Null verschieden, die b,(k=1,2,...,1) 
rein imaginar oder eventuell null und man erhdlt fiir (4.1) | 


Q(z) ay a, | a, | an | 


| | | . 
P(z) eben  e (4.3). 


Gemass dem Determinantensatz tiber die Naherungszahler und Naherungs- 


nenner eines Kettenbruches ist P(z) als n-ter Naherungsnenner von (4.3) gleich 
folgender Determinante: 


a+b) =z a, 0 0 0 0 0 
Pz) = 0 0 = b= ae 0 0 (4.4)5 
TEESE 10'S (0: 0)'e) 0! 0/18 ve hie. ‘eos SiBrze, 21 ie; e)ie™\s! (:[0!n) a) 0, 6: wish le) erteuiay emai e rattemete neon oe , 
0 0 0 0) Oia tn Ope he, ; 
0 0 0 O° Oe) ST aes ; 
5 


Aus (4.4) ist ersichtlich, dass P(z) das charakteristische Polynom des Eigen- 
wertproblems (N—z E) x = 0 ist, wobei die Matrix N folgende spezielle j 


' 
2 
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bische Form aufweist: 


| a, +2 Uae Gq 0 0 0 0 | 
oh ye beeen ub Onan 
Ns 0) Oo =-1 4, Gs OM (4.5) 
0 eee tas! 4) be ae | 
0 OP Spend = Na 
In der Matrix N sind die Gréssen a, (k = 1, 2,..., 2) reell und von Null ver- 
schieden, die b, (k = 1, 2,..., 2) rein imaginiir fee null. 


Ubertragung des Rittentsnts von Wall auf die Matrix N: Die Stabilitatsfrage 
kann nach (4.2) allein aus dem Realteil a, des ersten Elementes in der Haupt- 
diagonalen und den reellen Elementen a, (k = 2, 3, ..., m) in der oberen beglei- 
tenden Diagonalen beantwortet werden. Nach (4.2) gelten: 


Rom Gene whee S%e ath a Sing, we a'e 67s © & 0, ue 6 10! 0) mils We) 0) 0 ©) silt) 6.0) 6 0) @ 80.6) as fe) o6,/al ola 


"(Rk gerade) | 


(k ungerade) 
und es ist somit 


spnr, = —Ssgn (a; dy@3-+-a,) (k=1,2,...,n)- (4.6) 


Demnach folgt aus Satz 1: 


Satz 5: Die Zahl der positiven Glieder in der Folge von Produkten ay, a, 4p, 
Ay Ay Ag, .». , Ay Ay *** Ay _1 An ergibt die Anzahl der Eigenwerte von (N —z E) x =0 
mit positiven Realterlen. 


Schlussfolgerung: Das gestellte Problem, das charakteristische Polynom 
einer gegebenen Matrix A mit komplexen Elementen vermittels einer Kette von 
solchen Polynomen zu ermitteln, welche die Stabilitatsfrage beantworten kann, 
ist damit auf die Aufgabe zuriickgefiihrt, die Matrix A mit Hilfe einer Ahnlich- 
keitstransformation in die Normalform N iiberzufihren. 

Im zweiten Teil wird gezeigt werden, dass diese Transformation im allge- 
meinen durchfiihrbar ist. Dabei erweist es sich als zweckmiassig, die Normal- 
form N an der Nebendiagonalen zu spiegeln und die Bezeichnung leicht zu 
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Andern. Es seien 


Cb. Oe aan 0 0 
sob 4 Gee WUpge tne 0 0 
Be 0) ail bs Dea amen 0) 0) (4.7) 
0 0 0) 0 i Uysaa bya 
0) 0) 0) 0 ci. 1 by, n ar t Oa 


bzw. im reellen Fall 


Ot Vb TeeO 0 0 0 
ae Sees 0 0 
Be Oo Sa OM elas 0 0 (4.8) 

0 0 0 0 ees ace 

0 0 0) 0 —1 “ban 
die Normalformen der transformierten, zu A Aquivalenten Matrix, worin die — 
by, pay (R= 1, 2,...,—1) und 6,,, reell und von Null verschieden sind, wah- 
renddem die b;, ; fie =1,2,..., 1) in (4.7) reell aber nicht von Null verschieden — 


sein miissen. An Stelle von ban 5 lasst sich dann folgender Satz formulieren: 


Satz 6 (Stabilitdtskriterium fiir Matrizen): Wenn es gelingt, eine gegebene — 


Matrix A = (a;,) mit komplexen Elementen a;;, vermittels einer Ahnlichkeits- 
transformation in die Form (4.7), bzw. (4.8) bet reeller Matrix A, tiberzufiihren, 
dann ist die Zahl der positiven Glieder in der Folge von Produkten by n, Onn On—im 
Ban Onin On—-2n-1 +++ > Onn Onin °** O03 542 gleich der Anzahl der Eigenwerte von 
A mit positiven Realteilen. Sind speziell alle Elemente byy, bys, .. 
positiv, so haben alle Eigenwerte negative Realtetle. 


II. Methode der Transformation 


5. Elemente der Transformation 


Die Grundidee der Transformation besteht darin, durch geeignete Wahl und 
Reihenfolge von méglichst einfachen Ahnlichkeitstransformationen in endlich 
vielen Schritten sukzessive die gewiinschte Normalform zu gewinnen, indem i in 
jedem einzelnen Schritt ein Element der Matrix behandelt wird. 

Einerseits bewirkt ja Multiplikation mit einer Diagonalmatrix T von links 
TA in A eine Multiplikation der Zeilen mit den entsprechenden Diagonalele- 


of Oa came -Uin : 
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menten von 7, und Multiplikation von rechts AT’ eine entsprechende Multiph- 
kation der einzelnen Kolonnen von A. Anderseits bewirkt ja Multiplikation 
mit einer nicht Diagonalgestalt aufweisenden Matrix T von links TA in A eine 
lineare Zeilenkombination und Multiplikation von rechts AT eine lineare Kom- 
bination der Kolonnen. Als einfachste Gestalt wird T fiir einen Schritt einer 
Ahnlichkeitstransformation TAT-1 entweder als Diagonalmatrix so gewahlt 
werden, dass nur je eine Zeile und eine Kolonne in A verindert wird, oder als 
Nichtdiagonalmatrix so, dass nur zwei bestimmte Zeilen und Kolonnen mit- 
einander kombiniert werden. 


Typus A: T ist Diagonalmatrix 


Es sei 
1 eee 0) 0 Lives Oa 0) 
OSS « Oey O 0) Ome Oia 0) 
i Qomkfa 0 0 Ono Somer ) 
T = || Pe os 
ORO 1 Ol’ 000) 4 0 |? (Sai) 
OO (0 1 . Oumr O Ou) i 


wo 1/a in T, bzw. « in der Inversen 7-1, in der 7-ten Zeile und 7-ten Kolonne 
an Stelle von 1 steht. Eine Ahnlichkeitstransformation TA 7-1 mit einer solchen 
Matrix werde im folgenden mit 


| (Z 5%) | 


bezeichnet und hat auf A folgende Wirkung: 
Division der j-ten Zeile durch « und Multiplikation der j-ten Kolonne mit «. 


Typus B: T ist Nichtdiagonalmatrix 


Es sei 
1 ro 0 bet GG 0 
LO Wee 0 1 Ove 0) 
eae enor sok OS bee i OOO Dee WO? ete ie (Sez) 
Cae ere) 1 O-OFFOee 0 i 


wo ain T, bzw. —a« in T-}, in der j-ten Zeile und k-ten Kolonne steht. Die 
Matrix T, von links multipliziert TA, addiert in A zur j-ten Zeile das «-fache 
der k-ten Zeile und die Matrix T-1, von rechts multipliziert A 7-1, subtrahiert 
in A von der k-ten Kolonne das «-fache der j-ten Kolonne. Demnach bedeutet 
eine Ahnlichkeitstransformation TA T-1 mit einer Matrix T von diesem Typus: 
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Man addiere zundchst zur j-ten Zeile von A das «-fache der k-ten Zetle und | 
subtrahiere dann von der k-ten Spalte das a-fache der j-ten Spalte, oder auch in | 
umgekehrter Rethenfolge. } 

Eine solche Umformung kann als « zweidimensionale » Transformation be- 
zeichnet werden und sei im folgenden durch 


| (Lie: | 


dargestellt, indem die Indizes die Position von « in der dazugehorigen Matrix T 
und damit die zu aindernde Zeile (7) und Kolonne (k) sowie ihre Kombination 
kennzeichnen. 

Versucht man nun, die Transformation auf die Normalform mit Hilfe dieser 
beiden Typen von Umformungen durchzufiihren, so erkennt man bald, dass sie 
im wesentlichen in zwei Hauptabschnitte zerfallt: Im ersten Teil werden die 
Elemente unterhalb der Diagonalen behandelt, wahrend im zweiten Teil die 
Elemente oberhalb und in der Diagonalen auf die gewtinschte Form transfor- 
miert werden. 


6. Umformung der Matrix unterhalb der Hauptdiagonalen 


In diesem Abschnitt sei die Aufgabe gestellt, die gegebene Matrix A = (a;,) 
vermittels einer Ahnlichkeitstransformation auf die nachstehende Form zu 
bringen : 


Pu Piz Pis Pras a Pin 
mn Poo Pos Pos Sones Pon 
=e Og A Pss Psa Shea, p n 
Px 3 
0 Qos Pas Pan 6. 
‘ pi icaeas Corea 6 rae 


Dieses Problem ist allerdings bereits gelést durch das Verfahren von HEs- | 
SENBERG [8]. Doch beschreibe ich hier zur Vereinheitlichung der ganzen Trans- 
formation einen andern Weg. : 

Bezeichnungsweise: Zur Vereinfachung werden die Elemente der umge- 
formten Matrix immer wieder mit a;, bezeichnet, das heisst, es sei a;, der Wert 
des Elementes in der j-ten Zeile und -ten Kolonne, wie er sich nach der 
beendeten letzten Umformung ergeben hat. 


Die Transformation der Matrix unterhalb der Diagonalen kann gemiass 
nachstehender Vorschrift geschehen: 
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Regel 1: Beginnend mit der ersten Kolonne und von Kolonne zu Kolonne 
bis zur zweitletzten fortfahrend, fiihre man fiir die k-te Spalte (k= 1,2,...,n—1) 
folgende Umformungen aus: a nd 
| a) Um in der (k + 1)-ten Zeile eine —1 zu erhalten, wende man die Transfor- 
mation (T4143 —Ax41, 4%) an; 
| b) die Elemente a;, (j = k+2,k+3,...,%) werden dann mit Hilfe der zwet- 
dimensionalen Fak Jateibionen (Tia at zum Verschwinden gebracht. 
: Dass die Regel 1 den gewiinschten Erfolg hat, zeige ich folgendermassen: 
Voraussetzung: Die Elemente, durch welche dividiert werden muss, seien 
von Null verschieden. 
Behaupiung: Die Matrix A = (a;,,) mit komplexen Elementen a;, lasst sich 
unter Anwendung der Regel 1 auf die Form (6.1) transformieren. 


Beweis 


Induktionsannahme: Die erstenm— 1 Kolonnen weisen bereits die gewiinschte 
Form auf, und es seien also 


eae = —1, 46-0. fir k= 1, 2,16, — 1; (= kh+ 2, kh +3, 0 
Induktionsschluss: Anwendung der Regel 1 fir k = m 
a) Nach Voraussetzung Sei 444, m # 0; die Transformation (Ty413 —@m41, m) 
kann ausgefiihrt werden. Danach wird zuerst die (m + 1)-te Zeile durch —@n.44, m 
dividiert, wodurch in den ersten (m— 1) Kolonnen nichts gedndert wird, da ja 
Wmir,%— 0 ist fiir k = 1, 2,..: , m —1; in der m-ten Kolonne wird @m44 m ZU —1 
und in den folgenden Kolonnen werden die Elemente a,,,,,;, mit R= m-+1, 
m+ 2,...,n verandert. Dann wird noch die (m+ 1)-te Spalte mit —a,44, m mul- 
tipliziert, das heisst die Elemente der nachstfolgenden Spalte. Die der m-ten 
yorangehenden Kolonnen werden durch diese Transformation nicht geandert. 
b) Die Elemente a;, (j= m+ 2,m-+ 3,..., nm) werden vermittels der 
Transformation (Tj, m4} 4m) ZU Null gemacht. Danach wird zuerst das a; ,,-fache 
der (m + 1)-ten Zeile zur j-ten Zeile addiert. In den (m —1) ersten Kolonnen 
wird dadurch wiederum nichts geandert, da nach Induktionsannahme a,;, = 0 
m= 1,2, . — 1) ist sowohl fiir j = m+ 1 als auch fiir 7 = m+ 2, m+ 3, 
ee , , eas in eit Fall 7 = (m — 1) + 2 ist. In der m-ten Kolonne ist jetzt 
Gmst.m = —1, weshalb a;,, in Null itbergeht und in den folgenden Kolonnen 
erden die Elemente a;, (k = m+1,m+ 2,...,) geandert. Dann wird noch 
von der (m+ 1)-ten Spalte das a;,,-fache Seas sie Spalte (7 =m + 2) subtra- 
u nert, das heisst, es wird nur die nachstfolgende Kolonne verandert. 
- Damit ist gezeigt, dass sich das Verfahren gemass Regel 1 auf die m-te 
fortsetzen lisst, ohne dabei die m — 1 ersten Kolonnen zu andern. 
 Induktionsverankerung: Samtliche Uberlegungen des Induktionsschlusses 


alten sinngemiss auch fiir m= 1. 


tt 
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7. Diskussion 


Die Transformation stésst dann auf Schwierigkeiten, falls nach beendigter 
Umformung der (k — 1) ersten Kolonnen das Element @;,,,;, Null oder auch 
nur sehr klein geworden ist, durch welches im nachsten Schritt die (A + 1)-te 
Zeile der Matrix dividiert werden soll. 


1. Fall: Die Transformation kann auf die k-te Kolonne ausgedehnt werden 


Es sei @yi1, x = 0, bzw. sehr klein, aber mindestens ein Element in der ’-ten 
Spalte a;, mit k+2<7 =” von Null verschieden, bzw. nicht sehr klein, 
Durch gleichzeitige Vertauschung der j-ten mit der (k + 1)-ten Zeile und der 
entsprechenden Kolonnen kann erreicht werden, dass nun in der k-ten Spalte 
und (k +1)-ten Zeile ein von Null verschiedenes Element steht, ohne dass 
dadurch das in den vorangehenden Kolonnen schon Erreichte zerstort worden 
wire. An Stelle dieser Vertauschung kann auch die im Endeffekt gleichwertige 
zweidimensionale Transformation (T;,,;,;; 1) eingeschalten werden, indem die 
j-te Zeile zur (k + 1)-ten Zeile addiert und dann die (k + 1)-te Spalte von der 
j-ten Spalte subtrahiert wird. Jedenfalls wird in diesem Fall die Regel 1 aut 
die k-te Kolonne anwendbar. 


2. Fall: Die Matrix zerfallt 


Neben dem Element a,,,,, verschwinden saémtliche Elemente a;, in de: 
k-ten Kolonne mit k + 2 <7 <n, bzw. werden sehr klein. Dann zeichnet sich 
eine Unterteilung der Matrix in Untermatrizen ab: ‘ 


Ay, yg eee Ay ei Sree + a j 
1 agg wes Ageia Aan! Gonzr +++ Gan ‘ 
0 = 3x1 Age: Spear azn 4 
| i : 
es 0 0 Lo agp: Gnas Gin (7.1 
wee nae ee eee cece ee eee enn ns coer ew eww se nsmesn as aeneanene canes Pimecnnecnensactusnesenan sosewscesesss====--9508 7 
0) 0) 0 0 Fett kt ceva aiin 
(@) : 
0 0 0 Frronti +++ Anean 
0 0 0 O22 


Dann lasst sich das Verfahren auf die quadratische Untermatrix rechts unter 
weiter anwenden, und man erhiilt so schliesslich eine zu A Aquivalente Matri 

welche nicht die Normalform (6.1) aufweist, sondern aus Kastchen der Forn 
(6.1) aufgebaut ist, welche sich langs der Hauptdiagonalen aufreihen. Ei 


3 
> 
+ 
is 
Ss) 
3 
3 
nahn a ae 
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solches Kastchen kann gegebenenfalls auch nur aus einem einzigen Element 
bestehen, welches dann gleich einem Eigenwert von A ist. Oberhalb der Kist- 
chen wird die Matrix im allgemeinen mit von Null verschiedenen, jedoch be- 
deutungslosen Elementen ausgefiillt sein. Die Untermatrizen von der Gestalt 
(6.1), die sich so langs der Diagonalen ergeben, brauchen aber nichts gemeinsam 
zu haben. Es kann sein, dass sie je voneinander verschiedene Eigenwerte be- 
sitzen, teilweise untereinander gleiche Eigenwerte aufweisen, oder dass alle 
Eigenwerte einer Untermatrix auch Eigenwerte einer andern sind. Zu jedem 
Fall lassen sich ohne weiteres Beispiele konstruieren. In gewissen Fallen ist 
diese Aufspaltung rein zufalliger Natur und kann durch geeignete Zeilen- und 
Kolonnenvertauschungen in der Ausgangsmatrix vermieden werden. In aiidern 
Fallen ist sie durch das Koordinatensystem bedingt und kann durch Ubergang 
zu einem anderen System umgangen werden. Ganz allgemein ist diese Erschei- 
nung immer dann vorhanden und nicht zum Verschwinden zu bringen, wenn 
die Matrix A mehrfache Eigenwerte besitzt, und fiir einen solchen der Rang- 
abfall grésser als 1 ist, mit anderen Worten, wenn das Minimalpolynom von 
A nicht mit dem charakteristischen Polynom iibereinstimmt. So ist es unmég- 
lich, eine Hermitesche Matrix mit mehrfachen Eigenwerten vermittels einer 
Ahnlichkeitstransformation auf die Form (6.1) zu bringen. Das Versagen des 
Verfahrens ist demnach in diesem Fall durch die Struktur der gegebenen 
Matrix bedingt, worauf aber hier nicht weiter eingegangen werden soll. 


8. Beendigung der Transformation fiir reelle Matrizen 


Voraussetzung: Die Matrix A sei auf die Gestalt (6.1) transformiert, und die 
Kettenbruchentwicklung (2.1) fiir das charakteristische Polynom existiert. 

Dann lisst sich die Matrix noch mit endlich vielen Umformungen vom 
Typus B auf die Normalform (4.8) bringen. 

Geht man an die Lésung dieses Problems, so’ist man versucht, die ge- 
wiinschte Form einfach dadurch zu erreichen, dass man die notwendigen Ele- 
mente durch geeignete Transformationen Schiitt fiir Schritt zu Null macht, 
indem man von Kolonne zu Kolonne fortschreitet. Fir drei- und vierreihige 
Matrizen gelangt man auf diese Art auch wirklich ohne grosse Schwierigkeiten 
zum Ziel. Bei hoherreihigen Matrizen scheint dies vorerst auch zu funktionieren, 
und man gelangt einmal zu einer dquivalenten Matrix, welche zwar in den 
(n —1) ersten Kolonnen die gewiinschte Gestalt aufweist, jedoch in der letzten 
Spalte im allgemeinen nur von Null verschiedene Elemente enthalt. Will man 
diese noch zum Verschwinden bringen, so wird das schon Erreichte teilweise 
wieder zerstért. Aus diesem Grund ist die Reihenfolge der Umformungen ge- 
schickter zu wahlen, und so hat sich fiir reelle Matrizen eine Transformation 
herausgebildet, welche sich in drei Teile gliedert und die sich in die folgenden 


Regeln zusammenfassen lasst: 


Ces 


; 
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Regel 2: Beginnend in der ersten Kolonne und fortschreitend von Kolonne zw 
Kolonne bis zur zweitletzten, werden alle Elemente a,;, mit gerader I ndexsumme 
j+tkud1isjs k<n—1 zu Null gemacht, indem im allgemeinen Schritt 
die zweidimensionale Transformation (Tj, 413 4x) ausgefiihrt wird. 


Regel 3: Die Elemente a;, in der letzten Spalte mit gerader Indexsumme 
j+nund j <n—2 werden zum Verschwinden gebracht, indem die Transfor- 
mationen (T;,n; —4jnl4nn) angewendet werden. 


Hat man samtliche Umformungen gemiass den Regeln 2 und 3 durchgefiihrt, 
so liegt nun eine Matrix vor, in welcher in und oberhalb der Hauptdiagonalen 
mit Ausnahme von a, », alle Elemente mit gerader Indexsumme Null sind: 


0 ays 0 ay 0 Gig 5 os Sees 
—1 0 Ags 0 Ags a SiGe 8 
Wh sal 0 Asa 0 Gigs eureka 
Aw 0) Oa. 0 Ags Oth aa aa (8.1) 
0 0 0. ¢— 0 Ait. aseeeeeatee 
0 0 0 0 0 0 Dan 


Die restlichen Elemente werden noch wie folgt behandelt: 


Regel 4: Beginnend in der letzten Kolonne und dann nach vorn fortfahrend 
bis zur vierten Kolonne, werden séimtliche Elemente a;;, mit ungerader Indexsumme 
j+kund4<57+3<5k<n in Null iibergefiihrt, indem im allgemeinen Schritt 
dte Transformation (Tj, 4; —45%/4x4, x) ausgetibt wird. 


Anmerkung: Die in den Regeln 2 bis 4 angegebene Reihenfolge ist nicht die 
einzig mégliche, sondern sie kann in bestimmten Grenzen noch variiert werden. | 
So ist im Prinzip jede Folge zuldssig, die so beschaffen ist, dass in ihr jeder 
Schritt das bereits Erreichte unverandert lasst. Im Hinblick auf die Verifi-- 
zierung des Verfahrens wie auch auf die Programmierung fiir Rechenautomaten — 
erschien mir die gewahlte Reihenfolge als zweckmiassig. 


' 
' 


9. Beweis des Verfahrens 


V oraussetzungen 


ee 


I. Die gegebene Matrix A lasst sich vermittels einer Ahnlichkeitstransfor- 
mation auf die Gestalt (6.1) bringen. 


II. Die Kettenbruchentwicklung des Quotienten Q(z)/P(z) (2.1) fiir das’ 
charakteristische Polynom P(z) = |zE—A| existiert. 
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Behauptung: Dann lasst sich A durch eine Ahniichkeitstransformation in die 
Normalform (4.8) iiberfiithren. 
Beweis: Nach Voraussetzung I kann angenommen werden, dass A schon 
die Form (6.1) aufweise, und es bleibt zu zeigen, dass die Anwendung der 
Regeln 2, 3 und 4 zum Ziel fiihrt. 


1. Regel 2 


Induktionsannahme: Die ersten (m—1) Kolonnen scien bereits nach Regel 2 
behandelt worden und man habe oberhalb der Diagonalen: 


@,=0 fir 7+kgerade (1Sj7Skim—1Sn-—2). (9.1) 
Induktionsschluss: Um das Element a;,, mit 7 + m gerade und j = m zum 
Verschwinden zu bringen, werde nach Regel 2 zuerst das a;,,-fache der (m + 1)- 
ten Zeile zur 7-ten Zeile addiert. Da in der (m+ 1)-ten Zeile die Elemente 
4n41,% = 9 sind fir 1 =k =m —1, so wird dadurch in den (m — 1) ersten 
Kolonnen nichts geandert; da @,.4,m = —1 ist, wird a,,, in Null iibergefiihrt, 
und es werden in der 7-ten Zeile nur die Elemente a;, mit m+1<k <n ver- 
andert. Dann muss noch von der (m+ 1)-ten Spalte das a;,,-fache der j-ten 
Spalte subtrahiert werden. Da infolge der Form (6.1) in der j-ten Kolonne nur 
die Elemente a;;+ 0 sind mit z<7+1<m+1, werden dadurch nur die 


Elemente a,,,,,,1n der nachstfolgenden (m + 1)-ten Spalte in und oberhalb der 


Hauptdiagonalen betroffen. Insgesamt werden also nur Elemente in und ober- 
halb der Diagonalen verandert, die in einem spateren Schritt behandelt werden, 


-somit ist die Transformationsmethode auf die m-te Kolonne fortsetzbar. 


Zugleich ist auch ersichtlich, dass mit Ausnahme von @;,, alle andern Ele- 
mente in der m-ten Kolonne ungedndert bleiben, weshalb die Transformationen 


fiir die m-te Spalte unabhangig voneinander in beliebiger Reihenfolge durch- 
_gefiihrt werden kénnen. Ferner sind die Umformungen nach Regel 2 immer 
-ausfiihrbar, da dazu nur Multiplikationen und Additionen notwendig sind. 


2 


Induktionsverankerung: Samtliche Uberlegungen des Induktionsschlusses 
gelten sinngemass auch fiir 7 = m = 1. 


2. Regel 3 
Voraussetzung: Nach Anwendung von Regel 2 hat man 
a;,=0-fir 7+ & gerade ts. h5"7— bs (9.2) 


Hier ist Division durch a,,, erforderlich. Dies ist jetzt gleich der Spur der 


f gegebenen Matrix A und stimmt somit mit b,,,, in (4.8) iiberein. Dieses ist aber 


i 
i: 


4 
- 


- 
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als Folge von Voraussetzung II und Satz 1 von Null verschieden, die Division 
also durchfiihrbar. 

Um a;, mit 7 + gerade und 7 S” — 2 zum Verschwinden zu bringen, 
bestimme man zunachst den Wert aj, = —4jn/4nn und addiere zuerst das 
a;n-fache der n-ten Zeile zur j-ten Zeile. Die n-te Zeile enthalt nur die beiden 
von Null verschiedenen Elemente a, , in der letzten und —1 in der zweitletzten 
Kolonne. Demnach wird a;, zu Null, wahrend das Element aj, »4 (j < 7 — 1) 
verandert wird, welches aber von ungerader Indexsumme und deshalb im all- 
gemeinen auch noch von Null verschieden ist. Sodann muss noch das «, n-tache 
der j-ten Spalte von der n-ten Spalte subtrahiert werden. In der q-ten Spalte 
sind jedoch nach (9.2) nur Elemente a;; von ungerader Indexsumme und mit 
i <7 +1<n—1 von Null verschieden. Da einerseits 7 und m gleiche, ander- 
seits 7 und 7 aber ungleiche Paritat haben, so verandern sich in der letzten 
Kolonne nur Elemente a;,, mit ungerader Indexsumme oberhalb der Diagonalen. 

Damit ist gezeigt, dass unter Anwendung der Regel 3 die Elemente a;,, mit 
gerader Indexsumme und 7 < ” — 2 in Null tibergefiihrt werden, wobei sich 
insgesamt in der letzten und zweitletzten Kolonne nur die Elemente mit unge- 
rader Indexsumme oberhalb der Diagonalen andern. 


3. Regel 4 
Voraussetzung: Es ist 


4;,=0 fir 7+ gerade 1S) 5h 5245451). (9.3) 


Induktionsannahme: Man habe gemiss Regel 4 schon die Elemente a,, in 
den letzten (n — m) Kolonnen zum Verschwinden gebracht: 


4;,=0 fir 7+kungerade (15795 k-—3;m+15k Sn). (9.4) 


Induktionsschluss: Das Element a;,, in der m-ten Kolonne mit ungerader 
Indexsumme und 7 < m—3 soll gemidss Regel 4 behandelt werden. Dazu 
bilde man zunachst den Wert oj m = —@jm/@m_1, m unter der spater zu rechtfer-_ 
tigenden Annahme, dass a,,_,, » + 0 sei. Dann addiere man zuerst das %; m-fache 
der (m —1)-ten Zeile zur j-ten Zeile. Die (m —1)-te Zeile enthalt nach (9.4) 
nur zwei von Null verschiedene Elemente, namlich Am, m m der m-ten und —1 
in der (m — 2)-ten Kolonne. Somit wird 43m ZU Null, und a; »~». mit ungerader 
Indexsumme und 7 < m — 2 verdndert, welches in einer vorangehenden Ko- 
lonne oberhalb der Diagonalen steht. Dann muss noch das Oj m-tache der j-ten 
Spalte von der (m —1)-ten Spalte subtrahiert werden. In der j-ten Spalte sind 
nach (9.3) nur noch Elemente a, ; mit ungerader Indexsumme und mit 


15521 Saye 


oe” he 
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von Null verschieden. Da jetzt einerseits j und (m —1) gleiche, anderseits 7 und ; 
ungleiche Paritaét haben, so andern sich in der (m — 1)-ten Kolonne nur Ele- 
mente a;, ,,_, mit ungerader Indexsumme oberhalb der Diagonalen. 

Damit ist erwiesen, dass sich bei der Behandlung der m-ten Kolonne nur 
Elemente in der (m—1)-ten und (m= 2)-ten Kolonne mit ungerader Index- 
summe oberhalb der Diagonalen andern. Daraus ergibt sich, dass sich im spe- 
ziellen die Elemente a;,_, , mit m <k <n nicht mehr dndern, sondern bereits 
gleich den Elementen 6,_, , in der Normalform (4.8) sind. Diese sind nach 
Voraussetzung IT und Satz 1 von Null verschieden. Im besondern ist Cinig m0, 
womit die getroffene Annahme zutrifft. 

Induktionsverankerung: Simtliche Uberlegungen des Induktionsschlusses 
gelten sinngemiss auch fiir m = n. 


Zusatze 


1. Die Matrix zerfallt in Untermatrizen. Die Voraussetzung I ist nicht erfiillt. 
Wenn jedoch die Kettenbruchentwicklungen (2.1) der charakteristischen Poly- 
nome der einzelnen Untermatrizen langs der Diagonalen existieren, dann lassen 
sich die Regeln 2 bis 4 auf die einzelnen Kastchen anwenden. 

2. Die Kettenbruchentwicklung existiert nicht. Die Voraussetzung II ist damit 

-Aquivalent, dass die Sturmsche Kette fiir das charakteristische Polynom in der 

Form (3.1) mit linearen Quotienten (R, y + S,) existiert. Ist dies nicht der Fall, 
so zeigt sich dies dadurch an, dass entweder b,,,, oder ein Element b,_,,; in (4.8) 
verschwindet, so dass entweder nicht nach Regel 3 oder 4 weiter transformiert 
werden kann. Deshalb versagt unser Verfahren in denjenigen Fallen, in denen 
die Sturmsche Kette (3.1) nicht existiert. 


10. Verallgemeinerung des Verfahrens auf komplexe Matrizen 


Die Transformation der gegebenen komplexen Matrix A in die Gestalt 
(6.1) lasst sich genau gleich wie im reellen Fall nach Regel 1 durchfithren. Hin- 
gegen muss die Methode, aus der Form (6.1) die Normalform (4.7) zu gewinnen, 
dahin modifiziert werden, dass an Stelle der Matrix (8.1) zunachst eine solche: 

von der Gestalt (10.1) angestrebt wird. 


5 es CI OE rae NCE cari eer on arar tac cmc 
A — C oe 0 0 = i) 1 Cya C45 Fe RCM OM a i ee ee I (10.1) 
0 0 0 0 0 ee re eee 
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Darin sind mit Ausnahme des Elementes in der letzten Zeile und letzten 
Kolonne im oberen Dreieck die Elemente mit gerader Indexsumme rein ima- 
gindr oder aber reell, falls ihre Indexsumme ungerade ist. Infolge dieser Modifi- 
kation wird in jedem einzelnen Schritt entweder ein Realteil oder ein Imaginar- 
teil eines Elementes zum Verschwinden gebracht, und deshalb wird in der 
Transformationsmatrix T 


Lei: Ome tae sega |e O.crss pangl 0 

Oded Oomyy 0) [On A iaaesige 0 
Teen Guctene o}) Tt=||0.0 1 0 1. 0} ao.) 

Ge UO 1] O Neon aee 1 


y entweder reell gleich « oder rein imaginar gleich 7 8 gewahlt werden. Um dies 
zu unterscheiden, werde ich mit (T;/,,; «) die Ahnlichkeitstransformation TA i? 
mit reellem y =a und mit (T;/,; 8) diejenige mit rein imaginérem y =7 B 
bezeichnen, wobei der entsprechende Wert von y in T in der 7-ten Zeile und 
k-ten Kolonne steht. Weiter setze ich 

jy = Ap tt a5y, (10.3) 
und dies soll wieder der Wert des Elementes in der 7-ten Zeile und k-ten 
Kolonne bedeuten, wie er sich nach Beendigung der letzten Umformung er- 
geben hat. Dann lasst sich das Verfahren in die folgenden Regeln zusammen- 
fassen: 


Regel 2a: Beginnend in der ersten Kolonne und fortschreitend von Kolonne 
zu Kolonne bis zur zwettletzten, werden sowohl alle Realteile a;,;, mit gerader Index-— 
summe 1 +k und 1 575k <n—1 vermittels der Transformationen (Tj p43; 
%—=A5x) wie auch alle Imagindrteile aj, mit ungerader Indexsumme j + k und 
1sj<kSn—1 mit Hilfe von (Tj x44; B = a},) zum Verschwinden gebracht. 


’ 


Regel 3a: Die Realteile der letzten Kolonne aj, mit gerader Indexsumme 
7+ nund | <n — 2 werden zum Verschwinden gebracht, indem die Transforma- 


tionen (Tjn; &= —@j,/@y) angewendet werden. Die Imagindrteile der letzten j 
Kolonne. a;, mit ungerader Indexsumme j +n und 7 <n —1 werden in Null 
uibergefuhrt, indem die Transformationen (T;',; 8B = —a',|a',») zur Anwendung ; 
gelangen. ; 


: 
Regel 4a: Beginnend in der letzten Kolonne und dann nach vorn fortfahrend 

bis zur dritten Kolonne, werden sowohl sdmtliche Realteile ai, mit ungerader i 

Indexsumme 7+ k und 4<j7+3<k<n durch die Transformationen (Tj p43 

& = —A5,/A,_1,,) als auch alle Imaginérteile aj, mit gerader Indexsumme ee k 
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und 3<7+2Sk<n mit Hilfe der Transformaticnen (Farge Pies a falas a) 
in Null tibergefiihrt. 

Bei der Durchfithrung der Umforraungen hat man auch hier innerhalb der 
einzelnen Regeln bestimmte Freiheiten in der Reihenfolge der zweidimensio- 
nalen Transformationen. 

Dass die beschriebene Methode unter den gleichen Voraussetzungen wie im 
reellen Fall und unter strikter Anwendung der Regeln 1, 2a, 3a und 4a zur 
Normalform (4.7) fiihrt, lasst sich mit den genau gleichen Uberlegungen wie 
fiir reelle Matrizen beweisen. Der Unterschied besteht nur darin, dass -man 
erstens zwischen den beiden Arten von Transformationen (y reell oder rein 
imaginar) und zweitens zwischen reellen und rein imaginairen Elementen zu 
unterscheiden hat. Deshalb sehe ich von der Wiedergabe dieses evidenten 
Beweises ab. 


11, Erganzungen 


a) Berechnung des charakteristischen Polynoms P(A) 


P(A) =|AE—A|=|AE— B| lasst sich aus der reduzierten Matrix (4.7) 
bzw. (4.8) auf bekanntem Weg rekursiv wie folgt berechnen: Man beginne mit 


und bestimme dann die Polynome ;,(A) 
PrilA) = (A — 6 css, 241) PulA) + Oi, x11 Pea() (11.1) 


bzw. im reellen Fall 
Prsild) =A Px (A) aie Dy, k+1 Py-lA) Cee) 


fe=0,1,2,...,”—1). Dannist 


P(A) = Pal) + (-1)"7* Onn bna(A) - (11:3) 


b) Zur Wahl der Rethenfolge der Umformungen 


Die Freiheit in der Wahl der Reihenfolge der Umformungen nach Regel 2 
bzw. Regel 2a, kann numerisch dahin ausgenutzt werden, die Reihenfolge 
derart zu wahlen, dass die entstehenden Elemente moglichst kleinen Run- 
dungsfehlern unterworfen werden. So lasst es sich hie und da vermeiden, dass 
unvermittelt sehr grosse Werte entstehen. Numerische Versuche in dieser 
Hinsicht lassen vermuten, dass durch Anwendung der angegebenen Regel 2, 


bzw. 2a, die Zahlwerte der beeinflussbaren Elemente giinstig ausfallen, falls 


_ ae ee ST 
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man die Elemente pro Kolonne oben beginnend zu Null macht. Diese Reihen- 
folge weist jedenfalls den beachtlichen Vorteil auf, gegeniiber jeder andern An- 
ordnung der zweidimensionalen Transformationen am wenigsten Multiplika- 
tionen zu bendtigen. 


c) Zerlegung der Transformationsmatrix in Links- und Rechtsdrevecksmatrix — 


Die Ahnlichkeitstransformation der gegebenen Matrix A in die Matrix P 
von der Gestalt (6.1) geschieht nach Regel 1 in einzelnen Schritten, in denen 
die Transformationsmatrix entweder Diagonalform (k +122 in Regel 1a) 
aufweist oder eine besondere Linksdreiecksmatrix (j > k +1 = 2 in Regel 1b) 
ist. Da das Produkt von Linksdreiecksmatrizen und Diagonalmatrizen selbst 
wieder eine Linksdreiecksmatrix ist, so hat die Matrix L in der Transformation 


TAL eer (11.49 


folgende Form: 
Leg 0 0 0 
OF los, <0 0 0 
OF. Leg 0 
L= ae 11.5) 
O lyn lag Nay 0 c ) 
0 Dy 2 Ly 3 Ly 4 Ly n 


£ 


Anderseits erfordert die Ahnlichkeitstransformation von P in die Normal- 
form B gemiss den Regeln 2 bis 4 (2a bis 4a) lauter Umformungen mit Rechts- 
dreiecksmatrizen, deren Diagonalelemente gleich 1 sind. Somit ist die Matrix R | 
in der Transformation 


RP RR (11.6) 


eine Rechtsdreiecksmatrix mit 1 in der Diagonalen: 


We tymg arena Tin 4 
Oe Leena ares Von 

Eee je 1 754 ’3n (11.7) | 
Ow OAH ODO 1 


Damit zerfallt die Transformationsmatrix T in 


TAT =B (11.8) 
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auf nattirliche Weise nach (11.4) und (11.6) in ein Produkt einer Rechts- und 
einer speziellen Linksdreiecksmatrix: 


Teel, (11.9) 


d) Der Rechenaufwand des V erfahrens 


Der Rechenaufwand eines Verfahrens soll gemessen werden mit der Anzahl 
der notwendigen Multiplikationen und Divisionen. Um die Stabilitatsfrage 
einer reellen Matrix von der Ordnung zu beantworten, braucht es bekanntlich 
mit dem Verfahren von WEBER-VOETTER [7] und mit der Methode der Sturm- 
schen Kette [5] insgesamt eine Zahl von Multiplikationen und Divisionen, 
welche sich asymptotisch wie 2* verhalt. Die Abzdhlung der Operationen des 
hier beschriebenen Verfahrens zur Reduktion der gegebenen Matrix auf die 
Normalform (4.8) ergibt einen Rechenaufwand, der sich asymptotisch eben-: 
falls wie n* verhalt. In dieser Hinsicht sind also die beiden Wege gleichwertig. 


e) Zusammenhang des Verfahrens mit der Methode der Sturmschen Kette 


Ein enger Zusammenhang zwischen den beiden Methoden ist ersichtlich, 
wenn man auf die nachstehende begleitende Matrix (11.10) die Rechenvor- 
schriften anwendet: 


TO nD ee 
=t—o '6 ‘eer. 
C= | Oe A pares O.. <6 (11.10) 
0 ya) ORC, > 
0) 0) 0) —1 Cn—-1 


In der letzten Kolonne der Matrix C stehen im wesentlichen die Koeffizienten 
des charakteristischen Polynoms 


PU) SI* — 0,55 11 A ee + HDT GA (eG, 


und die Ausgangspolynome /j(y) und F,(y) fiir die Sturmsche Kette sind ge- 
_ mass (3.2) 

Bel Ve eg a Og PO U2, Ve, 
&(y) = = by a Guns gis — Cg yn—B 1 oe. 
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Nach Anwendung der Regel 3, welche die Elemente der letzten Kolonne 
mit geraden Indexsummen in Null itberfiihrt, erscheinen in der letzten Spalte 


an den Platzen mit ungeraden Indexsummen im wesentlichen die Koeffizienten 4 


des ersten Restpolynoms F,(y). In der zweitletzten Kolonne hingegen stehen 


nun an den Platzen mit ungeraden Indexsummen im wesentlichen die Koeffi- — 


zienten des Polynoms F,(y), welche aber durch c,,_, dividiert sind. Im itbrigen 
hat sich in der Matrix nichts geandert. 

Nun werde die Regel 4 auf die letzte Kolonne angewendet: Die bewirkt, 
dass in der zweitletzten Spalte bis auf einen Proportionalitatsfaktor die Koef- 
fizienten des Restpolynoms F;(y) in Erscheinung treten, in der drittletzten 


A hag Ea its 


Kolonne diejenigen von i(y), welche aber durch den Koeffizienten der héch- — 


sten Potenz dividiert sind. Fahrt man nach Regel 4 fort, so lasst sich allgemein 
zeigen, dass im wesentlichen und bis auf Proportionalitatsfaktoren fir k = n, 
n—1,...,5,4 in der k-ten Kolonne die Koeffizienten von F,_419(y) und in 
der vorangehenden (k — 1)-ten Spalte diejenigen von F,,_;,,(y) stehen. Fiir 
k = 3 dagegen steht in der dritten Spalte im wesentlichen der Koeffizient von 
F,,_,(y) und in der zweiten Kolonne derjenige von F,,(y). 

Somit stellen die Regeln 3 und 4 in diesem speziellen Fall eine besondere 
Form des euklidischen Algorithmus dar. 


12. Beispiele 


1. Reelle Matrix 


Unterhalb den Matrizen ist die angewendete zweidimensionale Umformung an- 
gegeben, welche zur nachsten 4quivalenten Matrix fiihrt. 


| WS teil Z 
Cal et een ahaa amie) 
gt eh eee 
2 On Su. 3 
Regel 7 
(Fy3 3) 
| 1 —9 —1 Z 
eet lh 1 0 — 0,333 333 0 | 
| —] —3 4 —1 
2 0 a 3 
(T39; —1) 
ile — 10,000 000 — 1,000 000 2,000 000 
pe Ree al — 0,333 333 — 0,333 333 0 
0) 1,333 333 4,333 333 — 1,000 000 
Z — 1,000 000 — 1,000 000 3,000 000 | 


(L493 2) 


sacieentnentniatiy ih intent miami tin diet ateaiadl dll 
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1  ~14,000000  ~1,000 000 2,000 000 || 
is |-1 ~ — 07333333 — 0,333 333 0 | 
| 9 3,333 333 4,333 333° —1,000000 | 
| O — 7,666667  __—41,666 667 3,000 000 | 


(F3; —3,333 333) 


| 1 ~14,000 000 3,333 333 2,000 000 
|| -1- — 0,333.333 1,111 110 0 
eT ee 4,333 333 0,300 000 
Lan® cee: F.G66867 5,555 556 3,000 000 | 
(T43; — 7,666 667) 
1 —14,000 000 18,666 667 2,000 000 || 
21" "0.333.333 1,111 110 0 | 
a ae 6,633 333 0,300 000 
= O 0 — 22,299 998 0,700 000 
(Ty; 22,299 998) 
1  —14,000 000 18,666 667 44,599 996 | 
|} -1 — — 0,333 333, 1,111 110 0 | 
| 4h a4 6,633 333 6,689 999 | 
| 0 0 a | 0,700 000 | 
Regel 2: 
(Ti25 i) 
Ou s , =754;333 333 19,777 777 44,599 996 
eit 0,666 667 1,111 110 0 
oe bs ee 6,633 333 6,689 999 
. 0 0 | 0,700 000 || 
(To; 0,666 667) 
|. Gh = = 4ass 333 29,333 337 44,599 996 
: | mat 0 5,533 334 4,460 002 
T thus Qe ee 2 7,300 000 6,689 999 
, |} 0 0 aa 0,700 000 
j 
(Ty4; 29,333 337) 
| | 0  —14,333 333 0 65,133 332 
—_ 0 5,533 334 33,793 339 
ee" Gokerse: ¥ 7,300 000 6,689 999 
Ke. 0 at 0,700 000 
4 (Tq; 7,300 000) 
@ —.= 14,333:333 0 65,133 332 | 
, | = 0 5,533 334 — 6,599 999 
4 Sains aE 0 11,799 999 
| 0 0 =i 8,000 000 
, Regel 3: 
{ ' (Tq; 0,825 000) 
i 
— 


4 


{ 
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0 — 14,333 333 () 76,958 332 
| =i 0 4,708 334 0 
= 0 aul 0 12,624 999 
| 0 0 =a 8,000 000 
Regel 4: 


(Tyg; — 6,095 710) 


QO 8237623 0 0 

| sy 0 — 1,387 376 0 
Bee 0 hea 0 12,624 999 
| 0 0 =A 8,000 000 


Da nun Dy, > 0, bg > 0, beg < 0, big < O ist, und nach Satz 6 in der Folge 
von Produkten by, > 0, Dag bag > 0, Dae Bag Dag < 0, Dag Daa 503 O12 > 0, drei positiv 
und eines negativ sind, so haben drei Eigenwerte von A positive und einer nega- 
tiven Realteil (vgl. dazu R. ZURMUHLI [8], Seiten 295-296). 

Entwicklung des charakteristischen Polynoms nach (11.2) und (11.3) ergibt: 


p44) =90, 
pune ts, 
pilA) =A, 
p,(A) = 2 — 8,237 623 , 
p,(A) = 489,624 999 A, 
b,(2) = A4+ 3,000 000 42 — 103,999 982, 
P(A) = p,(A) + (— 1)? 8,000 000 p3(A) , 
P(A) = A*—8 A343 224 76,999 992 2 103,999 982. 


Das exakte charakteristische Polynom lautet 
P(A) = A4—8 A343 424-77 1-104. 


Die Rundungsfehler in den Koeffizienten des linearen und konstanten Gliedes 


sind bei der verwendeten Rechengenauigkeit von 6 Stellen nach dem Komma 
zu erwarten. 


2. Komplexe Matrix 


— 0,0560 — 15,19002 0,2856+4,7820i1 —0,0207+ 3,4710% 0,0138 —1,7190z 

A= ||—9:0318— 6,45107 0,0095+0,4850i —0,0105+ 2,1100% —0,0153—0,3460% 
— 0,0428 + 2,4600% —0,1762—2,9980i —0,0859—17,6620i1 —0,4729—0,9820i | 

|| —0,0244-+4+ 1,1550¢ -—0,0617—1,3540i -—0,0391— 7,14904 — 0,3845 —1,05902 || 


Die Umformung nach Regel 1 bietet hier nichts Neues, weshalb diese Trans- 
formationen nicht wiedergegeben werden. Man erhalt: 


—0,0560—15,19007 —24,2869+2,1361i —4,0878+68,6431i 14,5697+ 4,96944 ||_ 


ait 4 0,0326—0,25741 7,5268+ 0,4708i  0,1335—0,4619: ||_ 
; = —0,3091—17,29461  0,4362—0,0571 || 
0 = —0,1844 — 0,6840i 


rt ig. gl he a 
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1. Regel 2a: 
(Tig; « = —0,0560) 
— 15,19007 — 24,2888 + 1,29997 — 4,5093 + 68,6167% 14,5622 + 4,99531 
ie -— il — 0,0234— 0,25747 7,5268+ 0,4708% 0,1355 —0,4619% 
0) - 1 — 0,3091 — 17,29467 0,4362 —0,05714 
0 0 -1 — 0,1844 — 0,68407 
(Tig; B = 1,2999) und (735; « = —0,0234) 
—15,19007 — 24,2888 — 2,3419+ 68,21507 14,6364 + 5,56232 
he — 1 — 0,25747 7,5340 + 2,16937 0,1233 — 0,46062 
0) - 1 — 0,3325 —17,29461 0,4362 — 0,05712 
0 0 -1 — 0,1844 — 0,68407 
(Tig; « = —2,3419); (Ty; B = 2,1639) und (Ti; « = —0,3325) 
—15,19007 — 24,2888 68,21502 15,0683 + 46,95971 
he -— 1 — 0,25747 7,5340 1,2115— 0,86072 
0 — 1 —17,29461 0,4975— 3,41017 
0 0 — 1 —0,5169— 0,68407 
2. Regel 3a: 
(74,; B = 90,8487); (73,3; « = 2)3419) und (74,; B = —6,5971) 
—15,19002 — 24,2888 — 22,6337% — 1096,5394 
-— 1 — 0,25747 5,1903 1£23,22934 
a 0 —1 — 10,6974i 68,9012 
0 0 -— 1 — 0,5169—7,28117 
3. Regel 4a: 
(Ty53 a= 15,9147) und (T3,; B = =1,7885) 
—15,1900% — 40,2034 — 23,0395% 0 
— 1 D53114 — 0,7656 0 
i— 0 = 1 —12,4859% 68,9012 
0 0 -— 1 — 0,5169—7,28117 


(Ti2; B = —30,0941) 


14,9041: — 442,6538 0 0 | 
muy _ 28,56314 =" 0.7656 0 
Bs 0 ay) 1248591 68,9012 
0 0 ay — 05169 — 7,2811i 


Hier ist by, < 0, bgq > 0, bg3 << 0 und by, < 0, und damit by <0, bas bag < 0, 
Dug Ogg Dag > 9, Dag Ose Dag O12 < 0, weshalb nach Satz 6 ein Eigenwert von A mit 
Beativent Realteil vorhanden ist, wahrend die drei tibrigen negativen Realteil 


aufweisen. 
Das charakteristische Polynom lasst sich nach (11.1) und (11.3) berechnen: 


P(A) =1, 
p,(A) = A—14,9041i, 
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pa(A) = (A+ 28,5631 4) p,(A) — 442,6538 Po(A) 

— J2413,6590i A— 16,9468 , 
p,(A) = (A-+12,4859%) Py(A) — 0,7656 p,(A) 

= 34 26,1449% 12—188,2569 A— 200,1856: , 
by(A) = (A+7,28114) pa(A) + 68,9012 po(A) 

— J4 +4 33,4260i 13— 309,7201 22 — 629,79191 1+ 289,9266 , 
P(A) = p4(A) +0,5169 p,(A) , 
P(A) = 44 (0,5169 + 33,42602) 43+ (— 309,7201 + 13,5143) 22, 

+ (— 97,3100 — 629,79191) A+ (289,9266 — 103,47597). 


Anmerkung: Das Beispiel wurde urspriinglich mit grésserer Stellenzahl durch- 
gerechnet, aber hier mit den auf 4 Stellen nach dem Komma gerundeten Werten 
angefithrt. 
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Summary 


Given the eigenvalue problem (A — A E) x = 0 for real or complex matrices 
A the number of eigenvalues A with positive real parts is determined without 
evaluating the caracteristical polynomial. A proceeding is developed here to 
transform the given matrix A into a reduced form by applying a finite series of 
elementary transformations upon the matrix. The elements of the reduced 
matrix allow immediately to solve the problem. 


(Eingegangen: 19. Juni 1956.) 
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Free Lateral Vibration of Beams of Variable Cross Section 


By Epwarb W. Suppicer, Princeton, N. J., USA), 
and Nazi J. Tarren, Beirut, Lebanon?) 


1. Introduction 


The problem of free lateral vibration of non-prismatical beams, i. e., the 
determination of the natural frequencies and the shapes of the deflection curve 
during vibration, is of current interest since gaps exist in the considerable 
work accomplished in the past in this field. KrRcHHOoFF [1]%) was the first to 
solve the governing differential equation for the case of the conical cantilever 
beam and for the wedge-shaped cantilever with small vertex angle. CAPELLEN [2] 
presented a series solution for the cantilever of triangular width; Morrow [3] 
utilized an integral method and gave several applications; NICHOLSON [4] 
studied the double cone and some solids of revolution; MononoBE [5] consid- 
ered the hollow truncated cone; and WRINCH [6] the solid circular cone (for 
the higher frequencies and modal shapes). Other contributions were made by 
Ono [7] and Hort [8]. The integral equation method was presented by SCHWE- 
RIN [9], TEMPLE and BIcKLey [10], [11], and WHITE [12]. 

This paper deals only with the approach by the differential equation of 
motion for beams with cross sectional areas and moments of inertia that.vary 
in an exponential manner. Solutions are obtained for beams with different end 
conditions, and for two specific examples. For simplicity only, the calculations 
are limited to the first and second modes. 


Notation 


x,y,z Rectangular coordinates, 

ie span of beam, 

A(x) cross sectional area, 

I(x) moment of inertia of cross section, 
Xe normal function for beam, 

xX the z-th normal function, 

i function of time, 

‘ b; circular frequency for the i-th mode, 


1) Princeton University. 
2) American University of Beirut. 
3) Numbers in brackets refer to References, page 519. 
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gravitational acceleration, 
time, 

weight per unit volume, 
modulus of elasticity, 


a 
is) . 


constants, 
the i-th root of an algebraic equation, 
radil. 


Pati ta 
& 
> 


= 


2. Formulation of the Problem 
The well-known [13] differential equation*) for free lateral vibration of a 


beam is a ary ne 2 
(Eset) eat ae ) 


conditions and the initial condition. For a beam, the boundary conditions are 
ends simply supported, fixed, or free. For the initial conditions we assume that 
y and y are known at t= 0. 
To find y(x, ¢) that satisfies equation (1) let 


y= X(x) T@); 


substitution of equation (2) into equation (1) gives ‘ 


Hy 
. 
. 
We seek solutions (x, ¢) that satisfy equation (1) as well as the 7 


A. OFF, 
wpe (EIS a) T+. MO 
or 
g 1 ee 
yA xX dx? 


RY eer 
es) a 3 


the first term is independent of ¢ and the second of x. For a solution of equa- 
tion (1), evidently each term in equation (3) is constant. Let 
WGP IE 
F gq = const = —p? (4) 
which gives 
T=Acospt+ Bsin pt (5) 


and from equations (3) and (4) we obtain 


ax oa 
av (BU gee ea (6 


a hae 
) The effects of rotatory inertia and shearing force are neglected. 
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Or 
d4X gi- ax . ay Oe. piv 
| oii Pt Oe hi EBS EO SQ ES 5 
i de a fapaceh (7) 


from which X(x) is to be determined, subject to the boundary conditions for 
any given beam. Since equation (7) does not allow an elementary solution in 
terms of A(x), J(x) and its derivatives, we shall assign certain expressions to 
A(x) and I(x) in order to simplify the solution. Equation (7) is the ‘eigenvalue 
equation’ and the problem is to determine the values of #? (the ‘eigenvalues’ or 
‘characteristic’ values) for which it has non-trivial solutions under the prescrib- 
ed boundary conditions. The corresponding solutions X,(«) are called ‘eigen- 
functions’. These solutions are not unique since if X; is a solution of equation (7), 
so is C X; where C is a constant. Since equation (1) is linear, we may write its 
general solution as 


i=1 
wherein the constants A; and B; are determined from the known initial con- 
ditions in the following manner. 
If X,,, and X,, are any two solutions of equation (6) we may write 
9 2 2 
qd? Seles ex aah 


Wake Sakcleben eee 


ES ES ia Bae: 
dx (! dx? at tap 


Multiply the first of these two equations by — X,, the second by X,,, add and 
integrate over x from 0 to L. Thus 


zr 


a (22, — 2) fA X_Xpdx= { [Xn 2 (i at) Xn ae (Ee) as| 
0 0 


ye ae ne 


0 0 
aa 
eT Ge | d | a dX my 
Am Ga (! dx* )|+ | ae Ge Van 
0 0 . 
£ 
a AX iy aX, =) 
=X, ae (7 ax? ) DRAG dx? 
0 0 i 
a7. d2X AX m a 
err (! a tan ae 
0 0 
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For a beam, each end (i. e., for x = 0 or « = L) is either simply supported, 
fixed, or free. At a simply supported end, the deflection and bending moment 
are zere, thus 

doar 


X=I5>=0, 


At a fixed end, the deflection and the slope are zero 


dX 
melee rortem aE 


At a free end, the bending moment and the shear are zero 


aX d a*xX 
ia ay {f a) 0. 


It follows that the terms on the right side of equation (9) all vanish. Hence 


L 
(b3,— B8) [A Xa X dae 
0 
If (p;, — ~2) + O then 
i 

[A Xm Xn dx = 0 (m+n). (10 
: 
To determine the A; and B, of equation (8) in terms of the given initial 
conditions, we assume that 


Equations (8) gives 


y= y X; (—f; A; sinp, t + p; B; cos; ?). weal 
4 
Thus ; 


2, 4eXi= le), SP Bee = ele). 


Multiply both these equations by X; and integrate 


Ub, 


L L L 
[Ar Xpax= [10 X,dx, [6B XPax = [9(x) Xj ax 
0 0 


0 


or 


: ; iB 
J ta) XxX, d% [a(x X, dx 
Sate 


. in 5 ul4 
0 
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i. e., the constants A; and B, can be calculated after an expression for X; has 
been found that satisfies equation (7) and the boundary conditions. Substitut- 


ing the A; and B; back into equation (8) we obtain the unique solution of 
equation (1). 


3. Solution for X 
To simplify the solution of equation (7) for X, we assume 
ey ST Ag ae a Ae oT (Xxx J, co = J, 0? (13) 


wherein 6 and ¢ are arbitrary constants (note that the radius of gyration of 
the cross section of the beam is constant). 
Substitution of equation (13) into equation (7) gives us an equation with 
constant coefficients 
atx ax 


aX pr yA 
65s 2 vi 
dae oo aa a dx 2 er, Soo (14) 


The general solution of equation (14) is 


A= y age 


4=1 


where 1; are the roots of the characteristic equation 


29 A 
r+ 2bro4 be — Ete - 0. (15) 
Let ; 
2 pe 16 
ak 4 Egl, (16) 
then equation (15) may be written as 
72+ br—K)(?+b7r+K)=0, 
A... 
r2tbr—-K=0 or #?+br+K=0 
and the roots are As 
—b+)Ve2?+4Kk 
oa . 
2 
et ya “le 
Sp=V4K 402 and S,=V4K- 0 (17) 
‘wherein 4 K > 62. The four roots are 
aha =i 5 2S Ss 
ee oe, ES = +; i= aa ae 2, tis mere 


-—_,. -~ ~~ & * 
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and the general solution of equation (14) may be written as 


‘ B23 2 
EN 275 #2 (a, e142 ae Ay e751 #12 of As e2 5242 =f aye di 


let 

a,=a+f8 and a,=a—f8 
then Sr 

ay eg #2 at As e751 42 te | (a a B) go 42 fe (a hs p) ea AEG 
f S 
=2acosh2* + 26 sinh 42 
and 
e524? — cos ond L ¢sin Shits Ss «Gig a = cee SE MP yn ; 

we may now write X as follows X= en beh 

: -S 3 en Se 
x (2 « cosh es + 2 6 sinh a + [a, + a4] cos 2s + [a3 — a4] 7 sin 5 ) 


Or 


X =e (C, cosh Sit Cysihh 12: C, cos oa eee *2*) (18) 
wherein C,, Cy, C3, C, are constants. 

In order to formulate the boundary conditions we shall need the following 
derivatives of X 


; dX Lacey a : O40 Se 
Xt = =e PSC, [S; sinh a — bcosh—4 | 


2 

Saw 
2 

Sa 


+ C, [Si cosh — bsinh a] 


+ C; [—Ss sin 


er [ Se cos aie b sin eae 


: 
| 
xX" a — goal? {C1 [- 6S, aah ve ft (S? + 5?) aes *34) 
é 
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a asx —b x/2 2 
= aie oy Ic, |- _ cosh — + —> sinh 72% 
2 Sig! ot aah b Sx 
Ae Pr? ‘ee 
+ (st + B) | F sinh- z ees cosh- 5 )| 
> Si 4/65 
j + Cy (3 + 6%) ( " cosh ee sinh >? =) 
fie 8 2 
} = >i sinh = -++ b 1 cosh Si =| 
21) 
wi Boe ( 
+ C, [=P cos 28% p> ain | 
S Saray See 
a ave: 2 2 
col "32) ( 3 Sin- 5 — +. = COs 37) 
+ Cy [- : 2 sin aa -} 2 cos Pat 
+ (b? — S3) (42 cos a ‘af 2 Lae =) | 


4. Applications 


We seek the frequency equation (the equation whose roots give us the 
frequencies of vibration of the normal modes) for beams with different end 


conditions. From these roots it is then possible to obtain the corresponding 
modal configurations. 


Case 1. Simply Supported Beam 


The boundary (or end) conditions are: 


9 x 
(X),.5= 9, (X),27 =0, N17. 7. 
Gs Mae] AS. tet ME ae 
Figure 1 


Simply supported beam. 
From equations (18) and (20) we thus get 


dae Oe Sal 5 ill 
cosh puss C, + sinh : Cy + os Cy + Sie Ce 04 


(22) 
a 
Siar. eteeeene® : Uh Rare ad ib 


f(L) Cy + g(L) Cz + A(L) Cz + 7(LZ) Ca = 0- 


& 
&. 

= 
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If all the C’ + 0 (i. e., the deflections + 0) then the determinant of the coeffi- | 
cients of the C/ in the above homogeneous set of equations must vanish. 
Expanding this determinant, we obtain 


Sill, Hee 
(Oise) ees sinh> — 8) DF Samos (cos a cosh ; 1) i= 0. = (Zag 


wherein, as before: S2?= 4K + b? and S3;=4K —6?. The unknown in the 
frequency equation (23) is K. The value of the circular frequency p follows then 
from equation (16). 

Equations (22) also give us the ratios 


G ] 

oes =] 

C, 
Si 5 ae _ Spl | 

C, DS (cosh 3 cos 3 | 4 K sin | 

~ = ie ’ ‘ 

C1 b (s, ani BERG sinh >) (24) 
Ls Sve : Side 

eS DSH (cos =8= — cosh ==) —4K sinh 

C b (Si sine = SSia EE 


=) ‘ 
Equation (23) is satisfied by an infinite number of values of K, i.e., an 


infinite number of p. Call these values £; and the associated quantities K;, S;,,_ 
See » B; , Cj1, Cia, C33 and C;,. The general solution of equation (1) becomes — 


y= 3 Co * (A; cosp,; ¢ + B; sin; #) | : 
(25) = 
: Shae C; Swot Cie eee 
x cosh ote Ag yall 43 PROM ie IE! g2 ) | 
( sph Ci gee Ge COS a eet Ca sin —7—}. 


The constants A; and B; are obtained from the initial conditions as given by 
equations (12), and may ‘aihsork the remaining C constant. 


Limiting Case: If b = c = 0 then from equations (13) we have A = A, and — 


I = I,, 1. ¢., the simply supported beam with constant cross section. Equation ‘ 
(23) reduces to: 


ne Sale = 
sin = 0 or sin) KES 6. 


the consecutive roots of this equation are 


Se en ee ee 
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Equations (22) give us: 

A = Cy wake 
m=O, C,=0, == as Coe 0 vand «C.ee 0, 


and equation (25) becomes: 


Le °) 


— > (A; cosp,; t + B; sin; t) sin! z omy 
j=1 
these results agree with the well-known solution [13]. 


We proceed next to solve the frequency equation (23) for 


t=) S.A = Age” and T= 1,2" -and*o= i 
Let 
k= KL? (26) 
and from equations (17) 


(S,=V4K += lV4k+1, S,=V4K—=-Vse-1. 


Equation (23) becomes 


+ 1)sin/e—2 sinJ/n > 
(8 k? — 1) sin|/k z sinh |/& PS 
16k 1 Wee ee ee 
+ V/16 k 1 (cos|/ 7 cosh |/t + Z 1) =|: 


Denote the left side of this equation | by f(k). Both parts of f() are positive, 


thus & is chosen so that sin Vk — (1 1/4) and cos Vk — (1/4) (1/4) have opposite signs. 
Trials show that k ~ 10, so that 


(27) 


8k2-1~8h? and V16k2?—-1—4hk, cosh |/ +4 Sk, ‘ete 
Equation (27) may be approximated by 
k tan /e- — = tanh e+ + — =-> (27a) 


and the graphical solution is indicated in Figure 2. 

The intersection points (1, 2,3,...) give us the approximate values of 
Vk — (1/4), i.e., the & values. 

For k = 9-77, equation (27) gives /(k) = +8-5. 

For k = 9-79, equation (27) gives /(k) = —13. We accept as the first root: 


ky = 9°78 = (3-12)2. 
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ile fl 
&TaNVR-7, fanh Laan 
Figure 2 
Graphical solution of frequency equation 


The second root is in the neighborhood of 27 in Figure 2, 1.e 


al 
= —1-68. 


For k = 39-4, equation (27) gives /(A) 
For k = 39-6, equation (27) gives /(k) = —0-1. We accept as the second root 


he = 39-6 = (6:29)2. 


From equations (16) and (26) 
72 Pi y Ao Ri. Egly 
K; EB gl, wwe OF P= je y Ay 
Thus the first two circular frequencies are 
(3-12)2 ee (620)? [Eee 
y Ay ) Po L2 y Ay ’ (28 


pi = L2 
it may be seen that 6-29 — 3.12 = 3:17 ~ a and Figure 2 shows that the root 


are spaced at intervals of approximately 7? for the higher modes. 
For the first mode (c = 1, 6 = 1/L, k, = 9-78, etc.) equations (24) vedio to: 


Ey C C. 
Sab BAG 3 wee Ed Der 
“is So ee 
and equation (25) gives us the first modal configuration 
Xe 
(29 


(cosh3 17 [— 1-115 sinh3-17 * — cos3-09 + + 7-43 sin 3-09 ) 


this is plotted in Figure 3. 
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For the second mode (c =1,b=I1/L, ky = 39.6, etc.) equations (24) reduce to: 


— = <i -_ = Cy — . 
Cy or r , C3 = 13 6 


and equation (25) gives us the second modal configuration 


{& 634). 4 63% 6:3 
ee iene (cosh _ ot  ginhs =e eae — $13-6 sin 23% ; (30) 
E, ue iD 
this is also plotted in Figure 3. 
Plot of equation(30) for x, | 
NandX, 2 A | x | X 
-Plot of tion (2g) 
_~ Plot of equation (29) for X, 0 0 0 
0-2 357 11-66 
x 0-4 9°56 7:22 
0-6 5:34 — 5-39 
0:8 3-19 — 8-79 
1-0 0-05 0 


Figure 3 
Modal configurations for Case 1. 


Case 2. Beam With Both Ends Fixed 


The boundary conditions are: a 
(i), = Deame eo. pods ae kala 


(2) be = 83 (X’) pO. Figure 4 


e Fixed end beam. 


From equations (18) and (19) we obtain 


oS. 


GATE 
, cosh ace sinh 23* Co + cos! ae oie sin Ces Ox 
ay (31) 


— bcosh 


= (-S, sin oe —b cos— 2) Cy + (Se cos 2 6 sin *s=) Cae ou. 
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Expanding the determinant of the coefficients of the C’s (compare Case 1), we 
obtain the frequency equation 


be ag 50 aes 
b? sin a2 sinh + S, Sp (1 cos 5 cosh 5 ) abe: (32) 


Equations (31) also give us the ratios 


Cis cos 2am cosh = | 

Cy. sinh ain 2 an ae ; | 

& a (33) 
tiene (cosh 5 Sep a | 

MS aah ale soa ae | 


The complete solution of Case 2 is equation (25) with the p; from equation (32) 
and the C-ratios from equations (33). 


Limiting Case: If b = c = 0, equations (13) reduce to A = A, and I= i, 
i.e. , the fixed-ended beam of constant cross section. Equation (32) becomes 


cosVK L cosh VK el 
as it should be [13]. 


As in Case 1, we proceed to solve the frequency equation for 
e=1 ».A=Ae and I=he* and b= 
and, as before 
= : oS eel il —_——_ —___ ——_—— EES 1 ——— 

k= KL. S,=V4K+ R= + yee+t, S,=V4K— b= 7 V4hoi. 
Equation (32) becomes 

eae re hat Ae Viera ras | eae es 
sin|/a— + sinh |/e + 1 + 16 2 — 1 i cos kk — = cosh b+ t)ao. 


(34) 
The first two roots are 


hy = 22-4 = (4-73)®, ky = 61-8 = (7-85)?. : 
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And the first two circular frequencies are 


_ (4:73)? )/EgT, M85)? |/E gl, 
p= |S py ES Ea (35) 
Note that 7-85 — 4-73 = 3:12 ~7, i.e., the adjoining ’s differ by about 22 
[as is evident in equation (34)]. 
For the first mode (c = 1, b = 1/L, k, = 22-4, etc.) equations (33) give 


kee io eon Sead, 
DORKS ag teers hse ek re 0-002 


and from equation (25) 


re x/2L (cosh4-75 + — 0-98 sinh 4:75 r — cos4-71 + + 0-99 sin4:71 4 


(36) 
this is plotted in Figure 5. 


For the second mode (c= 1, b= 1/L, k, = 61-8, etc.) equations (33) give 
i a ee Lae 
SS a hy ae ua ahs SC 


and from equation (25) 


X, = e 7 (cosh 7-88 + — 1-01 sinh7-88 _ — cos7-84 + + 1-02 sin 7-84 +] 


(37) 
this is also plotted in Figure 5. 


Plot of equation(36) for x, 


Figure 5 


Modal configurations for Case 2. 


Case 3. Beam Simply Supported at One End and Fixed at the Other End 
The boundary conditions are: 


=O, (X) Oe (Tx eee= 0; (Chen = 0. 


“=L 


(X) 


x=0 


ZAMP VII/32 


514 Epwarp W. SuppicerR and Nazin J. TALEB ZAMP 


a x 
VA 


~ n 
y 
Figure 6 
Beam with one end simply supported, other end fixed. 


From equations (18), (19) and (20) we obtain 


Geni 6.3 0 0) | 

cosh =e C, + sinh aa C, + cos 25 Cs 4 sie ey = 0; | 
LS Cbs Gye eee eee | (38) 
(s, sinh * — 5 cosh 25.) C, + (Sy eosh Ss! _ by sinh =) oo 

-- (S. sin ti + b cos 28) C3 + (S. cos oe b sin 28) Cy =e 


Expanding the determinant of the coefficients of the C’s (as in Cases 1 and 2), 
we obtain the frequency equation 


b Sy S, (1 — cos-2* cosh 3%) + 0% sin 72 sinh = 
4 4 = (39), 
ae (S, sin =e cosh ae Ss cos Sa? sinh "| =i 


Limiting Case: If b=c=0, the beam is prismatical and equation (39) 
becomes sum ee 
tan//K L =tanhVK L 
as it should be [13]. 


As in Cases 1 and 2, we proceed to solve the frequency equation for c = 1. 
Equation (39) becomes 


1642 —1(1—cos}/#— + cosn|/ + Spee ee r 
V cos |/ & — 7 cosh |/ & + 4} sin|/k q sinh|/.& + A 
+2k (Vae+r sin /e— 4 cosh|/# + 4 | (40) 
~V48—Teos\/a— 7 sinn|/e + 2) =o 

ai 4 = 


The first two roots are 


hy = 1655 = (4-07)2, ky = 51-0 = (7-15)?. 
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And the first two circular frequencies are 


(4-07)2 )/ Eel, 
fi = age or | 2 ps = - ad Beit Fgh 
y Ay ey y Ay 


Case 4. Cantilever Beam 


The boundary conditions are: 


Pes d rit x - 
(IX eno = OF ax (IX Vip = = 0, (X) v= 0, (X’),.. = 0. 
0 
= 
y 
Figure 7 


Cantilever beam. 


From equations (18), (19), (20), and (21) we obtain 
(Si + 0?) (b2—S 
2 2 


— b (3 S?+ b%) C, + S, (3 82 + S2)C, 


ly Sg oe LU gee eRe 0 


—— 


+ b (3 S2— 0%) Cy + S, (3 b?— S32 C,=0, 


Sik ¢,_0, 


= C, + cos hs 5 


C, + sin 


cosh 22 ~ C, + sinh na 


: TE, 
(S; Ue — cosh "3*) Cy + (S; cosh--2= — bsinh = 


— (S2 sin Aa b coat =) C3 + (Ss cos ade b sin 


Expanding the determinant of the coefficients of the C’s (as in Cases 1, 2, and 3), 
we obtain the frequency oe 


; VE 
‘oper (1 + cos Sa F cosh > ‘=| + 3 db? sin enh? 
2 2 2 (42) 
ee sewers . 
+25 (S; sin — cosh-> Behe S, cos 3 sinh ) == 0 | 


Limiting Case: If b = c = 0, the beam is prismatical and equation (42) be- 
comes < ty 
cosV/K L coshVK L = —l- 
as it should be [13}. 
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As in Cases 1, 2, and 3, we proceed to solve the frequency equation for c = 1. 
Equation (42) becomes 


es na ae oe aie 
16 Rb? = I (1 + cos|/ — = cosh |/k =f =) | 
Pa iia Sea Rai - 
+3 sin |/t a sinh /t + — 
4 4 | ( 
Vario onlfene 
+ 2{V4k + 1sin k p cosh ae 


+V4n—1cos|/e— 4 sinh /k + 1) =0, 


43) 


The first two roots are 
pea a (2:12) eee 2h eye 
And the first two circular frequencies are 


ie us | ee 


, (4-97)2 ee 
WAG : Es Ag 


5. Actual Sections of Surfaces of Revolution 


We consider now the actual shapes of beams that satisfy the assumed 


conditions 
Alg\ = Ave. TG) ie, 


(the square of the rad.us of gyration of the cross section of the beam). 
In general, the beams will be hollow, and appear as shown in Figure 8. 


Figure 8 
Beams with constant I/A. 


For hollow circular cross sections, we have 


AQ = Ay ON = (RO 7), T(x) = Tyee = A (Re a) Se 


ee fy ge eee R+% 
A Ay (0) 4 oe Pig 5” 
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Since 
ex/L 
R2 oes ye = a and R2 ail? v2 =< a ; 
s 41g 
We obtain 
2 J, Ay eo*lle aT A, ef til 
a= 0 0 é 
R? = ae “age ea hela Mie gear > ~ aes ; (44) 
But for 7=>0 
As log +” 40 
1 ey Az 
and #/L <1 
4n I, (R2 + 7) 
* os loa =a PASSO Oe 
= log aR or c Slog (R? — v2) ° (45) 


Equations (44) may be written as 


2 R® = (RQ +72) + (RR — 73) e°** and 272 = (R3+ 72) — (RR — 73) e°7/-. (46) 


We may assign any values to Ry and 7, provided we satisfy equation (45). 
Let: 
= Fo — const. 
Yo 


Equations (46) become 


. 


(Hagens gene and 2(tP Gen -@ er 
0 
(47) 
Equation (45) gives 
aa 
(ms log 5 = iy 
Or if we assign a value to c (say c = 1) then 
get 1) oe (@+ 1) 
et arty ro i! ae lea eat (48) 


Example 1 
For c= 1, equation (48) gives 7 < 2-16. Choose 7 = R, |v, = 1-46. Equa- 
tions (47) give 


a = 1:56 + 0:56" and (4) = 1-56 — 0-56 c=. 


ts 0 


Numerical values are given in Table I and the plot in Figure 9. 
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: 
Epwarp W. SupricER and Nazim J. TALEB wal 


ELLA LL LTT Zo 


CLLBLLLLLALLELELM LAL pf 


0-250 


C= 1, A=Ayge” r 


/L 


f 216 20, Saeco 


Ry Rg 2 p,; rad/s | f, c/s 


9-78 = (3-12)2 39-6 = (6-29)2... 


22-4 = (4-73)2 61-8 = (7-85)2... 


16:55 = (4-07)2 51-0 = (7-15)2... 


Bye 
J 2 
7 
4. 


4-5 = (2:12)2 24-7 = (4-97)2... 


Figure 9 
Summary and plot of Examples 1 and 2. 
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Table 1 
# 
7 0-1 0-3 0-5. 0-7 0-9 1-0 
R | / | 
1-479 1-524 1-576 | 2641 1-714 1:759 
0 
Y 
aa 0-971 0-897 0-799 0-660 0-427 0-197 
0 
R-y 
ee 0-508 0-627 0-777 0-981 1:287 1-562 
0 


Example 2 


Again c = 1 and we choose j = Ry/ry = 1-10. Equations (47) give 


eae = 1-105 + 0-105 e** and a + Pie oe 


0 0 


Numerical values are given in Table 2 and the plot in Figure 9. 


Table 2 

x 
— Ol 0:3 0-5 0-7 0-9 1-0 | 
L | | | | | | 
R 
= 1-106 1-118 kM GA esp so me 
Yo 
= 0-994 | 0-983 | 0-967 | 0-946 | 0-920 | 0-905 
‘9 

| Soe 0-112 | 0-135 0-163 0-200 | 0-250 | 0-276 
ra | 
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Zusammenfassung 


Es werden die freien Querschwingungen von Balken mit exponentiell ver- 
anderlichen Querschnitten und Tragheitsmomenten fiir verschied2ne Rand-_ 
bedingungen untersucht. Die Frequenzengleichungen, Eigenfrequenzen und 
Schwingungsformen werden angegeben (wobei die Auswertung einfachheitshalber 
auf die beiden ersten Eigenschwingungen beschrankt bleibt). 


(Received: April 15, 1956.) 
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Ein kompaktes Dekadenzahlgerat mit bedarfsweiser Umschaltung 
der Zahlkandle, Untersetzerstufen und Zahlwerke 


Von WERNER BeEuscu und DaNniEL MAEDER, Ziirich!) 


I. Anforderungen an Mehrfachzahlgerate 


In der Kernphysik stellt sich immer wieder das Problem, elektrische Impulse | 
aus Messapparaten separat, aber gleichzeitig zu zihlen: etwa die Gesamtzahl der _ 
im einem oder mehreren Strahlungsdetektoren erzeugten Impulse und gleic rT 
zeitig dazu verschiedene Zahlen von Impulsen, die nach gewissen Kriterien 
(Grosse, Gleichzeitigkeit mit anderen Ereignissen usw.) ausgewahlt sind. Anlagen, : 
die bis zu 20 verschiedene Zahlungen simultan durchftihren, werden sehr haufig 
gebraucht*). Bei Apparaturen dieser Gréssenordnung erfordert die Summierung 
mittels elektromechanischer Zahlwerke weniger Aufwand als eine rein elektro- 


1) Physikalisches Institut der ETH. 


2 7 ) 
) E. Barpincer und W. HAEBERLI, Impulsverstirker und Impulsspektrographen. Ergebn, 


exakt. Naturw. ’ 4 (1 53); A. B. VAN RENNES Pulse-A li y. 
La 
te 2 “ 3} ( ) mp itude Anal ae mM N Uc le ay Resew Ch, 
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nische Zahlung*). Die mechanischen Zahlwerke er! 


; me lauben eine rasche Ablesung 
insbesondere, wenn auf die Interpolation der Untersetzer verzichtet werden kann. 


Dies erfordert allerdings, dass man das Untersetzungsverhaltnis in jedem Kanal 
der jeweils vorhandenen Stosshaufigkeit anpasst, wobei im allgemeinen nur ein 
Bruchteil der eingebauten Untersetzerstufen ausgentitzt wird. 

Weiter ist es oft wiinschbar, zum Beispiel bei Messungen an radioaktiven 
Substanzen kurzer Halbwertszeit, ohne Unterbruch zu yihler aber mit Able- 
sungen in kurzen Zeitabstanden. i? 


Figur 1 


Frontansicht des umschaltbaren 5-bis-10-Kanal-Zahlgerats mit eingebautem Netzteil. 


Um eine gegebene Zahl von Untersetzern méglichst vollstandig auszuniitzen 
und um die ganze Anlage klein und handlich zu halten, haben wir dem Bau des 
in Figur 1 abgebildeten Mehrfachzahlgerats folgende Bedingungen zugrunde 
gelegt*): 

1. Verwendung von 10 dekadischen Untersetzereinheiten®). 

2. Jede Dekadenstufe kann wahlweise entweder 
a) einzeln betrieben, das heisst direkt mit einer Endstufe mit Zahlwerk ver- 

bunden werden, oder 
b) mit einer benachbarten Stufe zu einem hundertfachen Untersetzer kombiniert 

werden. 

Die fiir diese Umschaltungen beniitzten 5 Kippschalter sind in Figur 1 unter- 
halb der Dekadeneinheiten sichtbar. 

3, Jeder der 10 Eingiange fiihrt zunachst auf einen relativ hochauflésenden 
Zweifachuntersetzer, um die Zahlverluste in den nachfolgenden Dekadenstufen 


3) Etwa nach Art des von HuTCHINSON und Scarrottr beschriebenen Systems, welches bei 
erosser Anzahl von Kandlen aufwandmiassig giinstiger wird. Siehe Phil. Mag. 42, II, 792 (1951). 1 

4) Vorlaufige Mitteilung: W. Beuscu und D. MaeveER, Lin einfaches Dekaden-Zahlgerat mit 
umschaltbaren Kanal- und Untersetzungszahlen. Helv. Phys. Acta 29, Fasz. 3, 263 (1956). 

5) Der Firma Philips-Electronica danken wir fiir die anlasslich des ETH-Jubilaums erfolgte 
Schenkung von betriebsfertigen Steckeinheiten des Typs 88929/02 mit E11-Dekadenrohren. 


oumedan 
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} 
klein zu halten. In Figur 1 sind oberhalb der Eingangs- und Erdungsbuchsen — 
die zu den Zweifachuntersetzern gehorigen Interpolationslampchen erkennbar, 
Sie dienen nur der Funktionskontrolle. Fiir die nach 2b) zusammengefassten 
Dekadenpaare ist jeweils nur der geradzahlige (rechts befindliche) Eingang” 
wirksam. Sty % ; 

4. Zur Erméglichung ununterbrochener Messungen in fiinf Kanalen werden 
die 10 Untersetzer und Zahlwerke in zwei Gruppen (I und II) unterteilt. Ein” 
gemeinsamer Schalter (in Figur 1 links aussen) erlaubt i 
a) alle Kanale dauernd (siehe Figur 2: Stellung «I + II»), ; 
b) abwechslungsweise die eine oder andere Gruppe in Betrieb zu setzen, so dass” 

immer eine Gruppe zwecks Ablesung und Nullstellung stillgelegt ist. ; 

Soweit dabei die Dekadeneinheiten einzeln betrieben werden (20fache Unter- 
setzung), so wird jeweils der ganze Zahlkanal stillgelegt, so dass auch die Dekaden 
abgelesen und auf 0 gestellt werden kénnen. Sind zwei Dekadeneinheiten zu 
einem 200fachen Untersetzer zusammengefasst, so wird nur von den zugehGérigen — 
Zahlwerken eines stillgelegt. Wegen der fehlenden Interpolationsmoglichkeit ist _ 
diese Betriebsart auf hohe Impulsfrequenzen beschrankt. ; 

Entsprechend der alternierenden Betriebsméglichkeit sind zwei Riickstell- 
tasten vorhanden (siehe Figur 1, links aussen). ; 

5. Die Eingainge haben auf kurze negative Impulse (Anstieg < 0,5 us, 
Abfall > 1 ws) von =10V Amplitude anzusprechen, wie sie in Apparaturen — 
mit mehreren Zahlausgangen normalerweise verfiigbar sind oder sich allenfalls” 
aus Zahlimpulsen anderer Form mittels eines RC-Diodengliedes leicht herstellen 
lassen. ; 

6. Die Endstufen haben Zahlwerke von 5800 2 und 12 Hy (Anzug 10 mA, 
Abfall 2 mA) bis zu 20mal pro Sekunde zu betatigen. Der hohe Widerstandswert 
erlaubt, an Stelle eines Thyratrons eine normale Doppeltriode E90CC als End- 
rohre zu verwenden, mit dem Vorteil der unabhangigen Wahl von Stromleitungs- _ 
und Erholungsdauer. 


II. Eingangsstufen und Umschaltungen , 


Einzelheiten der Schaltung, insbesondere die fiir die unter 2. und 4. bespro- 
chenen Kombinationsmoglichkeiten notwendigen Umschalter sind aus dem 
Schema (Figur 2) ersichtlich, das eines von den 5 identischen Paaren aus je zwei 
benachbarten Zahlkanadlen darstellt. Auf den Einbau besonderer Diskriminator- 
stufen kann hier gemass 5. verzichtet werden, da die zu zahlenden Impulse _ 
praktisch stets aus Apparaturen stammen, welche alle nétigen Diskriminatoren — 
bereits enthalten, und es ist deshalb mdglich, fiir die Zweifachuntersetzung ein- 
schliesslich der Erzeugung der fiir die Dekadeneinheiten erforderlichen positiven | 
Impulse spezieller Form (Amplitude 13,5 V, Anstieg > 20 V/us, Abfall < 2 V/us) 
mit ever Doppeltriode je Kanal auszukommen. Nach Ansprechen des Zweifach- 
untersetzers wird die Kopplungskapazitat von 220 pF iiber den Spannungsteiler_ 
im Anodenwiderstand der einen Halfte der E90CC um 13,5 V umgeladen. Ein 
negativer Ladungsstoss wird sogleich (< 1 us) durch eine Diode 0A71 gegen die 
Gleichspannung von + 141 V abgefiihrt; eine positive Ladung dagegen wird nur 
langsam (iiber 0,7 MQ) gegen 0 V abgeleitet, und der Entladungsvorgang wird 
durch die Diode erst aufgehalten, wenn das Potential wieder + 141 V erreicht) 
was nach 14 us der Fall ist. Der so gebildete 14 us lange positive Sagezahnimpuls 
wird tiber eine zweite Diode auf den Eingang (Kontakt Nr. 2) der zugehorigen 
Dekadenstufe gekoppelt. Die zweite Diode hat den Zweck, die bei einem erneuten 
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A eaten des Zweifachuntersetzers wahrend der Dauer des positiven Sagezahn- 
Epes entstehende steile negative Flanke von der Dekadenstufe fernzuhalten 
wobei die totale Verdrahtungskapazitat des Kontakts Nr. 2 (rund 25 pF) ein 


+ 270V -19V +270V 


----------}---4 


+ 62V 


Eingangt ia sno E> D 


I fe. 


Ruckstellung GH1A19V) 


Eingang I 


E90CC 88929/02 E90CC 
Figur 2 
Schaltschema von 2 benachbarten Zahlkandlen. 


Widerstandswerte in kQ bzw. MQ, Kapazitaitswerte in pF bzw. uF. Der Gruppenumschalter 
«1/(1+1I)/L1» wird fii alle Kanalpaare mit einer gemeinsamen Achse betitigt. 


rasches Abfallen des betreffenden Potentials verhindert. Ohne diese Massnahme 
wiirde das Zuriickkippen des Zweifachuntersetzers grosse zusadtzliche Zahlverluste 
der Dekadenstufe verursachen (vgl. Abschnitt ITT). 

Die fertig erhaltlichen Dekadeneinheiten (Typ 88929/02) sind in Figur 2 nur 
durch das Schema ihrer Steckeranschliisse dargestellt. Zur Serieschaltung zweier 
Dekaden wird durch den doppelpoligen Kippschalter «x 20/ x 200» der Aus- 
gang (Kontakt Nr. 1) der unteren Einheit direkt mit dem Eingang (Kontakt 
Nr. 2) der oberen Einheit verbunden. In diesem Falle dient die zweite Diode des 
oberen Kanals zur Isolation zwischen dem oberen Zweifachuntersetzer (der zur 
Gruppe I gehért) und den vereinigten Dekadenstufen, wobei vorausgesetzt 
werden muss, dass der Eingang I keine Impulse empfangt, solange der Gruppen- 
wahler auf II steht. Bei der Zusammenschaltung eines Dekadenpaares muss daher 
der dazugehérige «Eingang I» geerdet und standig der «Eingang II» beniitzt 
werden; die Umschaltung I/II bezieht sich in diesem Falle nur auf die Zahlwerke. 

Anderseits miissen bei ununterbrochen abwechselndem Betrieb der eanzelnen 
Dekaden die beiden Einginge I und II miteinander verbunden werden, worauf die 


doppelte Ausfiihrung der Eingangsbuchse II hinweist, 
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Ill. Totzeitverluste 


2 os Sat aie mks Meee 


Folgt auf einen idealen Zweifachuntersetzer ein mit einer definierten Totzeit Ti 
behafteter «langsamer» Zahler, so ist die Wahrscheinlichkeit, dass der letztere i 


bis zum nachsten Schritt jeweils 


OIE E Ady? 
1 Impuls verliert: re [- pe ( 7 | , 
: AT)4 HIE \ 
2 Impulse verliert: ae E a - 


usw., wobei 2 = mittlere Frequenz der statistisch verteilten Eingangsimpulse am 
Zweifachuntersetzer. 


Der gesamte relative Zahlverlust unserer Dekadenstufen (T = 27 us) betragt 
also, falls die in Abschnitt II erwahnte Méglichkeit zusatzlicher Zahlverluste 
behoben ist: 


c= (Ad) peat [1+ AT cube: (AT)? (1+ = Beeb . 


2 3 Bre 4p 8) 


Die Totzeitverluste des wirklichen Zweifachuntersetzers (t = 1,5 ms) liefern einen 
weiteren merklichen Beitrag 
0, & At 


an die Gesamtverluste. Durch die maximale Zahlfrequenz der Endstufen und den 


es ee ee 


nr 


groésstmoglichen Untersetzungsfaktor werden die zulassigen Eingangsfrequenzen ~ 


auf 2 < 4000/s begrenzt. Mit obigen Parametern berechnet man zum Beispiel 


bei A =e LOO.0 2000 3000 4000 5000 Eingangsimpulsen/s 
Op Oy = 0,29 0,44 0,76 1,14 dos % Gesamtverlust 


Experimentelle Bestimmungen der Totzeitverluste mit Impulsen aus einem 
Szintillationszahler ergaben einerseits beim Vergleich der fiinf 200fach-Unter- 
setzerkanadle untereinander eine ausgezeichnete Ubereinstimmung; anderseits lie- 
ferten Kontrollmessungen mit einem schnellen 2!%fach-Untersetzer von bekannter 
Totzeit eine Bestatigung der theoretischen Verluste innerhalb + 0,5% bis zu 
Eingangsfrequenzen von 4000/s. Uberbriicken der Seriedioden in den geradzah- 
ligen Kanalen (vgl. Abschnitt II) bewirkte zusdtzliche Verluste von 2... 3%. 


IV. Endstufen 


Die ebenfalls in Figur 2 dargestellte Univibratorschaltung antwortet auf die 
relativ kurzen (14 ws) positiven Ausgangsimpulse der Dekadeneinheiten mit je 
einem Impuls maximaler Stromleitung durch das Zahlwerk wahrend 20 ms 
Dauer. Aus den in Figur 3 nach Oszillogrammen wiedergegebenen Spannungs- 
verlaufen erkennt man, dass etwa 10 ms spater die von der Selbstinduktion des 
Zahlwerks herriihrende Schwingung véllig abgeklungen ist. Der zwischen Zahl- 
werk und 0,05-yF-Parallelkapazitat liegende 1-kQ-Widerstand dient nicht nur 
zur Dampfung dieser Schwingung, sondern vor allem zur Erzeugung eines nahezu 
senkrechten, in der Amplitude aber begrenzten Abfalls an der Anode Al (linke 
Halfte der Doppeltriode in Figur 2) zu Beginn des Ansprechens, damit das 
Gitter G2 schon nach 1 us vollstandig sperrt. Von den beiden mit G2 verbun- 
denen Kondensatoren dient der kleinere (820 pF) zum Ubertragen des senkrechten 


2 ee ee ee te 
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Teils der negativen Flanke, wahrend der grossere (9,01 wF) zusammen mit dem 
Seriewiderstand von 0,2 MQ und einem weiteren RC-Glied dafiir sorgt, dass G2 
keine unzulassig grossen Spannungsspitzen erhalt und doch mdéglichst lange stark 
negativ bleibt. Die Kippdauer von 20 ms wird im wesentlichen durch die 0,01-uF- 
Kapazitat und den Ableitwiderstand von 4,7 MQ bestimmt; das erwahnte zu- 
satzliche RC-Glied (0,05 uF und 0,2 MQ) wirkt im Sinne einer Verlangerung der 


Kippdauer derart, dass der Schnitt der G2-Kurve mit der Sperrspannungs- Hori- 
zontalen schliesslich steiler wird. 


Zahlwerk-Univibrator 
6kQ,12Hy / E90CC 
0 —~s a mess 
-20 satel G1 (la: Seren ee Oe 


10 20 30 40 
— Zeit (msec) 


Figur 3 
Impulsformen an der Zahlwerkendstufe. 
G Gitter, A Anode. 1 und 2 bedeuten im Schaltschema (Figur 2) das links bzw. rechts gezeichnete 
Triodensystem. 


Mit fabrikneuen Zahlwerken lassen diese Endstufen eine Zahlgeschwindigkeit 
von 25/s zu; bei langerem Gebrauch empfiehlt es sich jedoch, unterhalb von 
20 Zahlschritten pro Sekunde zu arbeiten. 


V. Netzteil 


ti i r Dekadeneinheiten wird eine Anodenspannung von unge- 
a op re a ena ere : die Zweifachuntersetzer arbeiten aber bereits ee 
der Halfte einwandfrei. Wir betreiben die letztern deshalb eisai il Se 
0 V und der halben Anodenspannung, bzw. zwischen der halben ie 2 . a 
Anodenspannung (je 7,5 mA). Die Endstufen nehmen im Ruhezustand je : Se 
bei maximaler Zahlgeschwindigkeit etwa 7 mA auf, so dass der Pee abe ee 
des gesamten Anodenspannungsteils zwischen 160... 190 mA schwanke d 
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Aus Griinden der Gewichtsersparnis wurde der Netzteil (Figur 4) ohne Anoden- 
spannungstransformator ausgefiihrt. Um allzu grosse Wechselspannungskompo- 
nenten von den Betriebsspannungen fernzuhalten, muss bet der Inbetriebsetzung 
auf passende Polung des Netzanschlusses geachtet werden, wozu die Signal- 
lampe L dient. Da bei direkter Verbindung zwischen einem Pol der gleichgerich- 
teten Spannung und Netzerde und bei Vermeidung von Transformatoren nur 
Einweggleichrichtung méglich ist, muissen die Glattungsmittel um so reichlicher 


-o9 15A £220C 250 : 
| 2x FXCD36-22 4 —e 35....70mA dst +270V 
D3 
Ss S —e102mA Untersetzerstuféen 
OAH |S Forel ren mre 
etzerde S m m 
: eo Rall eae 1 lee 
ov 
6A 
eam 25] 350V 
+H62V, 30K] 
2A awL| ‘tS = es oe 
Untersetzer 7 ’ +162 V 
,-2x100 5602 /2W 
350V"_ == —>2mA +M1V 
E25 RSV 
—+32mA 
6,3V | 5K/5W 
ve ea 
6.3V 25u| 350V 
= i + 32mA 
Untersetzer I + 27V 
+NV bes 
fyscx | [Jaa 
-—5+60mA 7 
0,5 0,5 —35..70mA\ 
1 =e Pe 
all 9 9am RickstellungT 9 974 Ruckstellung I : 


= 19V. 


Figur 4 


Schaltschema des Netzteils mit direkter Kopplung (Vermeidung eines Anodenspannungstransfor- 
mators), Polaritatsanzeiger, HF-Entst6rung, Heiztransformatoren und Riickstelltasten I und II. 


dimensioniert werden; der in Figur 4 angegebene Kunstgriff, die Glattungsdrossel : 
mit Abgriffen zu versehen, erlaubte jedoch eine betrachtliche Reduktion des Auf- 
wandes. Die nach der Drossel vorhandene Brummspannung von etwa 5 V (von 
Spitze zu Spitze) spielt fiir den Betrieb der Zahlwerkendstufen keine Rolle; fiir 
die Untersetzerschaltungen wird die Betriebsspannung jedoch einem Drossel- 
abgriff (rund 15% der Gesamtwindungszahl) entnommen, welcher einen giinstigen 
Kompromiss zwischen zusatzlicher Schwachung der 50-Hz-Komponente (unter-_ 
kompensiert) und Zunahme der Oberwellen (iiberkompensiert) liefert. Die hoch- 
frequente Siebung mittels zweier Ferroxcube- Spulen D,, D, hat sich als hinreichend 
wirksam erwiesen, um das Auslésen der Untersetzer durch die Netzst6rungen 
vollig zu unterbinden. Wegen der betrachtlichen Warmeentwicklung der insge- 
samt 40 Rohren wurde das Gehause aus gelochtem Blech hergestellt, so dass sich 
eine kiinstliche Ventilation eriibrigt. ; 
Das ganze Gerat beansprucht wenig Raum, wiegt 16 kg und verbraucht bei 
maximaler Zahlfrequenz 160 W. Es wurde bisher zur Messung von Szintillations- 


er i ie wor 


rote, * 
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spektren eingesetzt und hat sich beim Verfolgen der Spektren kurzlebiger Akti- 
vitaten als besonders niitzlich erwiesen. Bei einiger Ubung lassen sich damit fort- 


heey Ablesungen in 5 Energiebereichen mit Abstanden von 10 zu 10 s durch- 
uhren. 


Summary 


A versatile pulse counter arrangement for nuclear radiation measurements is 
described. It consists of ten binary scalers (one E90CC tube per channel, dead- 
time 1-5 mus), ten commercial decade plug-in units (one E1T + E90CC, dead- 
time 27 ws), and ten mechanical registers (driven by one E90CC, 20 counts/s). 
The binary scalers accept negative pulses (= 10 V, rise < 0:5 us, decay > 1 ps) 
and serve as fast input scaling stages; they are designed to produce output 
pulses of the special shape required for triggering the subsequent decade stages. 

If all decade units are used singly, this equipment permits to count pulses 
in ten independent channels at a maximum rate of 400/s with negligible losses. 
On the other hand, one or several decade pairs can be made to work as scale-of- 
200’s, so that up to 4000 pulses/s/channel can be measured with counting losses 
< 1%. 

In cases where only 5 independent channels are needed, a switch common to 
all channels eliminates the dead-time otherwise required for reading out the 
counting results. This is done by alternative switching between even and odd 
registers, and in those channels which are connected as scale-of-20’s, also the 
corresponding decade units are stopped and can be read whilst the other group 
of channels is working. This feature has proved especially helpful in investigations 
on the scintillation spectra of short-lived radioisotopes. 

The power supply has been constructed without using an anode voltage trans- 
former, so that the whole apparatus is small and relatively light weighted. 


(Eingegangen: 22. Mai 1956.) 


Free Convection Past a Vertical Plate 


By Morton Fryston, Cambridge, Mass., WSA) 2) 


1. Recently there has been a noticeable increase in interest in free convection 
problems. Generally the methods used are one or another of the approximate 
approaches developed for classical boundary layer work. There are significant 
differences between the two however [1]%), and it may prove helpful to have 
exact solutions available. One such case is well known, i.e., the flow past a 
uniform temperature vertical plate where the equations can be reduced to 
ordinary differential equations. Often, however, an essential feature is that the 
plate temperature is not constant. It does not seem to have been noticed before 
that exact solutions can be obtained for a family of plate temperature variations. 
This note contains the derivation of the ordinary differential equations involved. 


1) Massachusetts Institute of Technology. ; 
2) This work was supported by the Naval Supersonic Laboratory, Department of Aeronautical 


Engineering, Massachusetts Institute of Technology, Professor Joun R. Markuam, Director. 
3) Numbers in brackets refer to References, page 529. 
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2. The free-convection flow past a vertical plate is governed by the following 


equations and conditions: 
Ub Ug +0 Uy = OU Uy + 8 P(E Ti) a 
Uy + vy = 0, 
uT,+vT,=aTy,, 
at y= 0) w= w=0, Dees, 
asy—->oo: u>0, T=T7,, aconstant. 


To put these equations into dimensionless form we set 


2\1/8 _ v2\1/8 _ , T 
u=(vg)au, x= (=) ae y= (=| y, v= (vg)8a, 0=—. 


§ & 


Then 
U Uy + U Uy = Uy, + BT, (8-1), 
Wiz +0,=0, 
and , 
a O54 0 0, == 0555 


where Py (= v/a) is the Prandtl number. 


Now it is well known [3] that these equations reduce to ordinary differential 


equations when the surface temperature is constant. Then we set 


y= oP ty aN), O= 1+ cg hly 2), 


ZAMP 


(10) 


where u = y,, V = —Yy,. It is natural to ask if there are any other surface tem- 
perature distributions which permit exact solutions. The most obvious generali- 


zation of equations (10) is to set 
y= A x? fly x9) , = 1S Boh ae 
(For convenience the bars are omitted over all quantities.) This leads to 
A* x2P+20-11(b + g) f'?— pf fi] = A Pt8Cf" 4 Bah, 
since B 7, = 1. Thus we must have 
4p 29-1L=fPt+3.4¢» O@=Pt3¢) oo p=1aag, Os Nee A ie 
Finally, with A=1, 


i 3 
2 —— ar fio SET), 


Cp eps oe ula 
o f Go) eee 
When i A aA a a 


i> 10o7, i Sor 


(11 


) 


(12) 


(14) 


(17 


These equations include the constant surface temperature example (o=.0 
as a special case. The equations can be solved for various combinations of o, B 


) 
) 
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and Py, although only a restricted range for w is of interest. Since q is usualiy 


negative, the leading edge is at the temperature of the medium far from the plate. 
Equation (11) now reads 


y= PAC +3)/4 tly x-V/4) : B= + Bx” h(y wp —1)/4) : (18) 


y= BHO NE Ary ylo-Did) | (19) 


and the heat transfer at the surface is proportional to. For instance, it is constant 
when » = 1/5 and n= yx ¥®. 


Finally, solutions to equations (15), (16), (17) can be used to test approximate 
methods and in stability calculations. 
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Zusammenfassung 


Diese kurze Mitteilung zeigt, dass fiir den Fall natitirlicher Konvektion eine 
Familie von exakten Lésungen gefunden werden kann, welche die regierenden 
Gleichungen zur Darstellung bringen. 


(Received: April 10, 1956.) 


Ein Elektronenvervielfacher als Pulslichterzeuger 


Von REINHOLD GERHARZ, Washington, D.C., USA?) 


Aufbau des Vervielfachers 


Als Breitbandverstarker haben Elektronenréhren mit Sekundarvervielfachung 
der Hochfrequenztechnik neue Anwendungsgebiete er6ffnet. Auch in vielen Zwei- 
gen der Impulstechnik hat diese Rohrenbauart Eingang gefunden, und der ein- 
stufige Réhrentyp EFP 60 wird in vielen Versuchsschaltungen benutzt. 

Um die Eigenschaften eines nach Angaben von Bay ee) entworfenen mehr- 
stufigen Vervielfachers mit Gittersteuerung zu untersuchen, wurden in diesem 
Laboratorium einige Versuchsréhren hergestellt. Ferner stellte die Firma « National 
Union» (Orange, N. J.) einige technisch verbesserte Rohren des gleichen Entwurfs 
zur Verfiigung, die als Typ HTIFA in den Handel gebracht werden. HosTETLER 
publizierte konstruktive Einzelheiten iiber diese Réhren, welche jedoch einige 


Erganzungen erfordern. 


1) George Washington University, Elektronisches Laboratorium. ae. ‘ 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 535. 
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Der Gesamtréhrenleitwert (G) betragt etwa 0,1 A/V und ist definiert durch 
den Quotienten von Stromanderung AJ,,, und Gitterspannungsanderung AU, 


Al, 
G= — 2, 
AU, 

Der Entwurf sah urspriinglich einen maximalen Anodendauerstrom (Kollektor- 
strom) von J,= 200 mA vor, der aber wegen der Warmebelastung des Elektroden- 
materials nur auf 100 mA gehalten werden konnte. Dabei war ein Betriebsfaktor 
(«duty-cycle») von 50% in Rechnung gestellt worden. ; 

An den héheren Dynodenstufen des Vervielfachers sind zur Erleichterung der 
Warmeabfuhr Abstrahlflachen angebracht, deren Oberflachen aufgerauht wurden. 
Als Dynodenmaterial wird eine Legierung von 98,3% Silber und 1,7% Magnesium | 
verwendet. Das Formieren einer sekundaremittierenden Schicht erfolgt durch 
Erhitzen der Dynoden in einer stark verdiinnten Sauerstoffatmosphare. Metal- ¢ 
lisches Magnesium (Mg) an der Dynodenoberflache wird dabei zu Magnesiumoxyd ~ 
(MgO) umgeformt, welches dann als diinne Schicht von einigen Angstrémeinheiten 
die Dynode bedeckt. Da das restliche, noch in der Dynode befindliche nichtoxy- 
dierte Mg bei hohen Betriebstemperaturen der Rohre an die Oberflache diffundiert 
und dort abdampfen kann, wurde als obere Grenze der zulassigen Dynodentem- 
peratur 500°C gewahlt. Auch bei dieser Temperatur ist die Mg-Verdampfung iiber 
gréssere Zeitabschnitte immer noch nachweisbar, und es empfiehlt sich, die Be- _ 
triebstemperatur unter 500°C zu halten. 

Der Hauptanteil der anfallenden Warmeenergie kann nur durch Strahlung aus 
der Rohre abgefiihrt werden. Deshalb lasst sich in grober Annaherung die Tem- 
peraturbilanz der Elektroden mit einer bei der Bestimmung der Konstanten o in 
der Stefan-Boltzmannschen Strahlungsformel angewandten Beziehung erhalten{[3]: 


4 W : 1 
Z Vee 
dabei bedeuten 
T Temperatur der sekundar emittierenden Oberflache, 
W Gesamtenergie des elektrischen Stromes durch die Dynodenoberflache, - 


o Stefan-Boltzmannsche Konstante, 
fq Verhaltnis der warmeabstrahlenden zur warmeerzeugenden Oberflache, 
€, relatives Emissionsvermégen infolge der Beschaffenheit der Oberflachen, 
&) relatives Emissionsvermégen des Oberflachenmaterials, 

I, Temperatur der Umgebung (hier des Rohrenglaskolbens). 


Es konnte festgestellt werden, dass eine wirmeabstrahlende Oberflache von 
tile 6 cm? geniigt, um die letzte Dynode unter der kritischen Temperatur zu 
alten. 

_ Der Ersatz des Silbers als Hauptmetall des Elektrodenmaterials durch eine 
Ni-Rh-Mg-Legierung?) unterdriickte die Verdampfung des Mg nur unwesentlich, 
wahrend die im folgenden zu beschreibenden Fluoreszenzeigenschaften der Dynode 
bei dieser Legierung merklich schwacher erscheinen. 

Die bei den hohen Stromdichten im Flugraum zwischen den Dynoden ent- 
stehenden storenden Raumladungseffekte, welche die Elektronenflugzeit erhéhen 


3) Ni als Trager wegen der besseren relati issi i 
) elativen Emission; darauf eine dii ic 
Verhiitung der Oxydation des Ni wihrend des F ‘ on altnatleceveaea eam 


anes ormierens. Beim Erhitzen diffundiert fast nur Mg 
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und die Sekundarausbeuten erheblich verringern, konnten durch Hilfselektroden 
fast ganz beseitigt werden. Dazu wurden bei diesem Réhrentyp Stabchen aus 


Tantal (@ © 0,8 mm) verwendet, die mittels Halteplattchen aus Glimmer unge- 


fahr in der Brennlinie der nahezu halbzylinderférmig gekriimmten Dynoden Auf- 


stellung fanden. Die Betriebsspannung fiir die Stabchen wurde von der tibernachst 
hoheren Elektrode abgezweigt. Im Falle der letzten Dynode liegt also die zuge- 
horige Hilfselektrode auf Kollektorpotential. Die durchschnittliche Sekundaraus- 
beute dieser Versuchsréhren konnte bei Stufenspannungen von 250 V und mit 
Hilfselektroden vor der 5., 6. und 7. Dynode auf 6 = 4 gehalten werden. Da die 
Raumladungseffekte nicht véllig zu unterdriicken sind, verursachen diese immer 
noch St6érungen in den Elektronenbahnen und dussern sich hauptsachlich durch 
Zeitverzogerungen von etwa 10-®s zwischen dem Eintreffen der ersten und der 
letzten Elektronen eines Signals auf der nachsthdheren Dynode. Daher kann die 
Zeitaufldsung eines mit dieser 7stufigen Réhre gebauten Hochfrequenzverstarkers 
nicht besser sein als 7 - 10~® s, was bei sinusférmigen Schwingungen einer Frequenz 
von 71,5 MHz entspricht. Die Gesamtflugzeit der Elektronen durch die Réhre 
betragt etwa 3,3 -10-*s und bewirkt bei hochfrequenten Schwingungen nur eine 
Phasenverschiebung. Die Eingangs- und Ausgangskapazitaten betragen bei dem 
vorliegenden Vervielfacher etwa 18 pF. 


Betrieb als Impulsgenerator 


Die dynamische Kennlinie des Vervielfachers fiir Eingangsimpulse von 1078s 
Dauer ist in Figur 1 wiedergegeben. Der Steuerteil ist in Aufbau und Wirkungs- 
weise einer Pentode vergleichbar. Bei der Kennlinie ist jedoch der mittlere Teil 
wegen der Elektronenvervielfachung um ein Mehrfaches steiler. Die weitgehende 
Nichtlinearitat der RGhrenkennlinie zwischen Anlaufstrom- und Raumladungs- 
gebiet wurde von Bay als giinstige Voraussetzung fiir die Erzeugung von rechteck- 
férmigen elektrischen Impulsen hoher Frequenz erkannt [1]. Um die aus der Rohre 
gelangenden verstarkten Impulse mit einem Zeitverzug wieder zum Steuersystem 
zuriickzubringen, wurde ein Koaxialkabel als Riickkopplungsglied verwendet 
(Figur 2). Die Lage des Réhrenarbeitspunktes ist so zu wahlen, dass von den an- 
kommenden Signalen nur solche ausreichender Amplitude in den steilen Teil der 
Kennlinie gelangen. 

Auf diese Weise lassen sich niedervoltige Rauschspannungen von dem erwtinsch- 
ten Impuls trennen. Die Ausgangsamplitude der Impulse ist infolge des horizontal 
auslaufenden Teils der Kennlinie im Raumladungsgebiet begrenzt, und die Impuls- 
dauery wird von der als Differenzierglied arbeitenden Kabelankopplung festgelegt. 

Bei einer Kabelriickfiihrung von der letzten Dynode zum Steuergitter laufen 
im Kabel positive Impulse um, wahrend bei der Riickfuhrung zwischen Kollektor 
und Kathode negative Impulse durch das Kabel laufen. In beiden Fallen betrug die 
Impulshdhe etwa 7 V und die Impulsdauer etwa 8 - 10-* s, was einer Impulslange 
auf dem Kabel von etwa 

Tam hp pits? Vinspilie = Pe S220 Sir 107% 200-0m 
entspricht. Die Impulsdauer lasst sich durch am Ende kurzgeschlossene Kabel- 
stiicke («Stubs» von L ~ 50cm) bis zu etwa 3° 10s verkiirzen (vgl. Fig. 2). 

Die Impulsfrequenz hangt von der Lange des Riickfiihrkabels ab und betrug 

etwa 10 MHz bei 25 m und 100 kHz bei 2500 m Kabellange. Das Riickftihrkabel 


a 
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ist zur Vermeidung schadlicher Reflektionen am Niederspannungsende mit seinem 


Wellenwiderstand abgeschlossen. 
J 
eae ea 
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Figur 1 Figur 2 
Dynamische Kennlinie eines 7stufigen Elek- Riickgekoppelter Impulsgenerator mit Reflek- 
tronenvervielfachers; Dauer des Eingangs- tionsstiick zur Impulsverkiirzung. Die Dy- 
impulses 10~8s, nodenstufenspannungen betragen 250 V. 


D 3 Hilfselektroden vor den letzten Dynoden; 
A Hilfsanode; C Kathode; H Heizung; L Riick- 
fiihrkabel als Laufzeitelement; U, Gittervor- 
spannung; X Induktiver Spannungsabfall und 
Kabelankopplung; Z Kabelabschlusswiderstand. 


Dynodenlumineszenz 


Beim Betrieb der Roéhre ist (besonders an den Dynoden mit héherer Strom- 
dichte) eine blauliche Lumineszenz der elektronenbeschossenen Oberflache zu 
beobachten. Die Intensitat dieser ins Sichtbare reichenden Strahlung scheint 
linear mit dem auf die Dynoden treffenden Elektronenstrom anzuwachsen und ist 
beim Impulsbetrieb der Leuchtstarke einer kleinen Kontrollglimmlampe fiir elek- 
trische Apparate vergleichbar. Die gesamte Lichtenergie ist zu etwa 0,2 W ab- 
geschatzt worden. 

_ In Hochvakuum-Vervielfacherréhren ist das Leuchten der Dynoden schon seit 
einiger Zeit bekannt. SCHAETTI [4] erwahnt bei seinen Photomultipliern mit hoher 
Verstarkung eine blaue Emission der Silber-Magnesiumplatten, welche aber durch 
die Verwendung von Cu-Be-Dynoden wesentlich verringert werden konnte. BENE- 
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= [5] berichtet von ahnlichen Leuchterscheinungen und versucht eine Erklarung 
urch die Annahme von Restgas-Ionisierung auf den Elektroden. Andere Autoren 
[6], [7], [8], [9], (1 0] untersuchten die Sekundaremission unter verdnderlichen phy- 
sikalischen Bedingungen und bringen Exkiarungsversuche; jedoch sind Leuchter- 
scheinungen wahrscheinlich wegen der geringen Stromdichte in ihren Apparaturen 
nicht beobachtet worden. Eine Anzahl von Arbeiten iiber geziichtete Leuchteffekte 
fiir technische Zwecke (Elektrophosphoreszenz) mit verhaltnismassig langen Nach- 
leuchtzeiten (10-6s) wurde kiirzlich wahrend eines Kolloquiums in England 
diskutiert [11] (vgl. auch [12)). 

Die in der untersuchten Vervielfacherréhre auftretenden Leuchteffekte wurden 
mit einem Spektrographen photographiert (Platte vom III F-Typ). Selbst bei 
einer stundenlangen Exposition konnte keine Banden- oder Linienstruktur des 
Spektrums nachgewiesen werden. Die Auflésung betrug bei diesem Gerat 15 A/znm. 
Die Wellenlangenverteilung der Lumineszenz erstreckt sich, aus dem Ultraviolett 
kommend, von etwa 3300 A an mit wachsender Intensitat iiber 4360 A (hier 
sanftes Maximum) und erreicht bei 4900 A schnell abfallend die Intensitat Null. 

Um die Dauer des Leuchtens zu untersuchen, wurde dieses mit einem 1 P21- 
Photomultiplier und mit einem Oszillographen (Typ Tektronix 517) bei einer Zeit- 
auflésung von 5-10-°s beobachtet. An das Gitter des Vervielfachers wurden 
mit einem Relaispulser (Spencer-Kennedy-Typ, etwa 100 Impulse pro Sekunde) 
positive Impulse von 7 - 10-*s Dauer gelegt. Dies entspricht der Zeitauflésung des 
Vervielfachers. Mehrere Messreihen zeigten, dass die Lumineszenz ohne beobacht- 
baren Zeitverzug den elektrischen Impulsen an den Dynoden folgt und dass die 
Abklingzeit des Pulslichtes nicht langer als 10-8 s dauert. 

Bei kleiner werdendem Abstand von Lichtquelle zum Photomultiplier konnte 
die statistische Amplitudenverteilung der vom Photomultiplier abgegebenen 
Impulse beobachtet werden. Vergrésserung der einfallenden Lichtintensitat bzw. 
Abstandsverringerung liessen erkennen, dass die einzelnen Impulse statistisch 
immer haufiger in einer Enveloppe von der gleichen Aniplitudenfunktion erscheinen, 
welche den elektrischen Impuls am Ausgang der Vervielfacherréhre beschreibt. 

Da die Lumineszenz nur in der diinnen MgO-Schicht auf der elektronenaufneh- 
menden Seite der Dynode beobachtet werden kann, wird die Abklingzeit des 
MgO-Leuchtens mit der gemessenen Abklingzeit identisch sein. Elektrolumineszenz 
von messbarer Intensitat, die auf reinen Metalloberflachen auftreten k6énnte, ist 
wegen der sehr geringen Lichtabsorptionstiefen der Metalle unwahrscheinlich. 

Infolge der kurzen Anregungs- und Abklingzeiten und wegen der emittierten 
Wellenlangen scheint es also berechtigt, das beobachtete Leuchten mit ‘Elektro- 
fluoreszenz im Innern des MgO zu erklaren. Als Voraussetzung fiir die Sekundar- 
elektronenemission der MgO-Schicht wird die dort herrschende hohe elektrische 
Feldstarke angesehen, welche je nach der MgO-Schichtdicke bis zu 1 08V/cm betragen 
kann. Jacops [7] erklart die hohe Sekundarausbeute durch lokale Bildung von 
Elektronenlawinen in den MgO-Kristallchen. Es scheint deshalb sehr gut méglich, 
dass ein Teil der so gebildeten Elektronen als Leuchtelektronen in freien Energie- 

bandern Austauschwirkungen hervorrufen und so die Fluoreszenz verursachen. 
. Die verschiedentlich vermutete Ionisierung von Gasresten an der Dynoden- 
flache als mégliche Ursache der Lumineszenz [5] ist wegen deren kurzen Abkling- 
zeit und der zu erwartenden, aber nicht nachweisbaren Struktur des Spektrums 
‘abzulehnen. Auch eine spurenhafte Reduktion des MgO und dabei freiwerdender 
Sauerstoff infolge der Elektronenbeschiessung kann wegen der hohen chemischen 
Affinitat der Komponenten zueinander nicht als Hauptursache des Leuchtens an- 
_ gesehen werden. Kurzzeitige Uberbelastung der Dynoden ergaben zwar Erhitzung 
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bis zur Rotglut und einsetzende Mg-Abdampfung, jedoch die Réhreneigenschaften 
erschienen nach dem Erkalten unverandert. 

In den meisten gegenwartigen technischen Anwendungen der Elektronenver- 
vielfachung werden Lumineszenzerscheinungen als nichterwiinschte Nebeneffekte 
behandelt. Der beschriebene Vervielfacher erlaubte es aber, gerade die Leucht- 
eigenschaften fiir mehrere Versuchszwecke auszunutzen: Eine Intensitaétsmodula- if 
tion des Dynodenleuchtens konnte zum Beispiel bis zu einer Frequenz von 50 MHz 


eae gee 


nachgewiesen werden, wenn sinusformige Hochfrequenzspannungen an das Steuer- 
system der Réhre gelegt wurden. Die Verwendungsméglichkeit der Elektrofluo- 
reszenz zu stroboskopischen Untersuchungen ist naheliegend, und auch die 
Erzeugung von Réntgenstrahlimpulsen hoher Frequenz und genau festlegbarer 
Modulation und Steuerung ist prinzipiell méglich. 


Pulslichterzeuger mit optischer Ruckkopplung 


Im Verlauf der Untersuchungen des Dynodenleuchtens erschien es aussichts- 
voll, bei dem in Figur 2 gezeigten riickgekoppelten Impulsgenerator das Riick- 
fiihrkabel teilweise durch einen Lichtweg zu ersetzen. Dies ist in Figur 3 schematisch i 
dargestellt. Alle Schaltelemente sind die gleichen wie in der Figur 2..Die vom ~ 
Vervielfacher emittierten Lichtimpulse wurden mit einem 1 P21-Photomultiplier F 
aufgenommen und iiber einen Verstarker von 10° Hz Bandbreite in das Riickfithr- 
kabel eingeleitet. 


Figur 3 

Oberer Teil des Vervielfachers mit Hilfselektroden, optischer Riickfiihrung und Laufzeitkabel L. ; 

M Verschiebbarer Spiegel; PM Photomultiplier; S Sammellinse ; D, Leuchtdynode; K Kathoden- 
verbindung. 


{ 


Der Lichtweg wurde mit einem Spiegel verandert. Eine Verschiebung des 
Spiegels von + 1cm entspricht bei dem gegebenen Kabel von 50m und einer 
urspriinglichen Impulsfrequenz von 5 MHz einer F requenzanderung von + 840 Hz. 
Ein Lichtweg von etwa 1 m bei 30 m Kabellange markierte bei der beschriebenen 
Schaltung die obere Frequenzgrenze (etwa 8,5 MHz), an der noch liickenlose Puls- 
lichterzeugung méglich war. Mit einer Sammellinse von 10 cm Offnung und 20 cm 
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Die Form der erzeugten Lichtimpulse muss in erster Naherung mit den elek- 
trischen Impulsen im Kabel tibereinstimmen. Dies wurde experimentell bestatigt 
Wegen Streulicht, Stérrauschen und Spannungsschwankungen gelangen in den 
Riickfiihrkreis eine betrachtliche Anzahl von Stérsignalen mit kleiner Amplitude 
Die impulsformenden Eigenschaften des Vervielfachers sieben jedoch auch in 


dieser Anordnung die Storsignale wegen deren begrenzten Amplitude im Anlauf- 
stromgebiet der Kennlinie aus. 


Diskussion 


Die Dauer der Lichtimpulse hangt ausser von der Zeitauflosung des Verviel- 
fachers nur von der Abklingdauer des lichtemittierenden Materials ab. Im Falle 
des MgO betrug die gemessene Lichtpulsdauer nur etwa 10-8s, und es konnten 
«Lichtpakete» erzeugt werden, die einer raumlichen Lange von 


L=cr=3-10- 10-8 = 300cm 
entsprechen. 

Sorgfaltiger dimensionierte Réhren und héhere Betriebsspannungen wiirden 
die Zeitaufl6sung um 2 Gréssenordnungen verbessern kénnen, wie es bei einer 
Variation des Vervielfacherprinzips von STERNGLASS und WACHTEL [13] kiirzlich 
angedeutet wurde. Hierbei handelt es sich um einen Vervielfacher mit mehreren 
planparallel iibereinander angeordneten Emissionsschichten (Dicke etwa 1000 A, 
mechanisch gehalten von einem Traggitter, Abstande voneinander etwa 1 cm, 
Stufenspannungen 2 bis 3 kV), in welchen die ausgelésten Sekundarelektronen 
durch die diinnen Vervielfacherschichten hindurch zur nachsthéheren Dynode 
gelangen. Die errechnete Zeitaufldsung soll 5-10~11s betragen. 

Abgesehen von stroboskopischen Anwendungen diirfte die Erzeugung von 
schnellen Lichtimpulsen auch fiir den Bau von automatischen Rechenanlagen 
bedeutungsvoll werden. In dieser Richtung arbeitend, hat LoEBNER [14] tiber 
Untersuchungen an «opto-elektronischen» Schaltungen berichtet, bei welchen 
elektrophosphoreszente Substanzen mit lichtempfindlichen Zellen gekoppelt waren 
und mit denen einfache Verstarker- und bistabile Systeme bei geringstem Raum- 
bedarf demonstriert werden konnten. Die -begrenzte Ansprechgeschwindigkeit 
dieser Anlagen im Millisekundenbereich behaftet praktische Ausfiihrungsformen 
jedoch mit schwer zu iiberwindenden Nachteilen. Bei Hochgeschwindigkeits- 
Rechenanlagen ist eine Loslésung der Informationselemente von materiegebun- 
denen Fiihrungssystemen in wachsendem Mafe erforderlich. Streckenmodulierte 
Lichtpakete eréffnen neue Méglichkeiten, ein Maximum von Informationselementen 
in eine gegebene minimale Zeitspanne einzufiigen und durch Regeneration fiir 
beliebig lange Zeit zu speichern. Ein Pulslichtverfahren kénnte die Anwendung der 
endlichen Ausbreitungsgeschwindigkeit elektrischer Wellen im freien Raum fiir 
die Speicherung von Informationselementen in die unmittelbare Zukunft riicken*). 
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Summary | 
7 


A 7-stage electron multiplier with grid control system is described with respect . 
to construction and performance. A coaxial feedback cable permits the regeneration _ 
of pulses with a repetition frequency of from 10 MC to 0-1 MC and a pulse length ~ 
of about 8 x10-®s. The luminescence of the dynodes was investigated and 2a 
identity in the duration of the light pulses and the electrical pulses has been found. — 

The interpretation of the luminescence effects by assuming electro-fluorescence . 
phenomena indicates applications to problems of stroboscopy and of high-speed 
computation and storage devices. i 


fom ae 


—, 


(Eingegangen: 29. Marz 1956.) 


Die Nusseltsche Zahl fiir laminare Str6émung im zylindrischen Rohr 
mit konstanter Wandtemperatur 


Von SERGEJ PaHor und JANEZ STRNAD, Ljubljana?) 


Fir die W&armeleitung in strémenden Fliissigkeiten gilt die Differential- 
gleichung?) 


OT 
y?T= io (5 v grad T). 


Fur stationare Strémung im zylindrischen Rohr (Durchmesser 2 7) mit konstan- 
ter Wandtemperatur Ty suchen wir eine Lésung der Form 
T= Ty + C R(x) exp(—£ 2) , (1) 


wobei #% = yr/v, und C = const. Damit erhadlt man fiir Poiseuillesche Str6mung 
die Gleichung 


R"(x) 4 =i) = ( _ ii + p (1 — x)| R(x) = 0 (2) 


me der Randbedingung R(1) = 0. Als Parameter erscheinen hier die Pecletsche 
ahl Pe= 207%, 9¢,/A und p= Bro Pe. Alle Stoffwerte wurden als konstant 


1) Universitat. 


2) H. GroBER, S. Erx, U. Gricuri, Grund 


; ; esetze der Wdrmeii i - 
Berlin, Géttingen, Heidelberg 1955), S. 179. i SS ae 
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vorausgesetzt. Hauptsachlich interessiert uns der kleiuste Eigenwert ~, und die 
zugehorige Eigenfunktion R,(x), welche in jedem Fall das Temperaturprofil fiir 
grosse Entfernungen von der Eintrittsstelle darstellt. 

In technisch wichtigen Fallen ist Pe su gross, dass das Glied (p/Pe)?, welches 
der Warmeleitung in der Langsrichtung entspricht, vernachlassigt werden darf. 
Die so vereinfachte Gleichung lésten Grartz und Nu SSELT?) mit Reihenentwick- 
lung. Hier wurden fp, und R,(x) mit iterierter numerischer Integration nochmals 
bestimmt. Es ist p, = 7,3140. 

Es erscheint uns mathematisch interessant, das Problem auch fiir beliebige 
Pecletsche Zahlen zu untersuchen. Fiir Pe +0 geht (2) in die Besselsche Dif- 


36 a 1 Penne | See Ss Ee 


"0 ‘) 10 15 20 


Pe=2V 1% 0C,/A. ——— 
Die Abhangigkeit der Nusseltschen Zahl von der Pecletschen Zahl. 


ferentialgleichung iiber, so dass R,(¥) = J)(k, ¥) und p,/Pe = k, = 2,4048. Mit 
k,, bezeichnen wir die m-te Nullstelle der Funktion Jy. ; 
Von beiden extremen Fallen ausgehend, erhalt man durch Perturbation 


Or125 
pi = 7,3140 (1 For pe -) (Pe > 1); (3) 
Py = 2,4048 Pe (1 — 0,1626 Pe + ---) (4a) 


und 
R,(«) = Jo(Ri *) + Pe [0,026 Jo(k, x) — 0,002 Jr(Rg *%) t+] + +++ (Pe <1). (4b) 

Die Nusseltsche Zahl wird definiert*) mit 
2g 27, (0T/0r) w 


Nu = 7 = Ty — TF 


, 


wobei 


Yo 
if 
. Te) = sag [Te?) v(r) 2avdr. 
0 


3) H. Groser, S. Erk, U. GRIGULL, Grundgesetze der Warmeiibertragung (Springer-Verlag, 
‘Berlin, Géttingen, Heidelberg 1955), S. 182. 
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Fiir ein Temperaturprofil der Form (1) ist 


1 
Nu = 73 


Mit den aufgestellten Ausdriicken (3) und (4) fiir p, und R(x) ergibt dies 


WBS, 
Nu = 3,6570 (1 dis Pe? =e -) (Pe > 1) 3 


und 


Nu = 4,1805 (1 — 0,0439 Pe+---) (Pe<1). 


Einige Werte von Nu fiir mittlere Pecletsche Zahlen wurden durch Polynom-_ 
approximationen fiir R,(¥) bestimmt. Den Gesamtverlauf der Funktion Nu(Pe) 
stellt das Diagramm auf Seite 537 dar. 


Summary 


The Nusselt modulus for laminar flow through a pipe with constant tempera- 
ture of the walls has been determined as a function of the Peclet modulus. 
This function is represented by a diagram. 


(Eingegangen: 11. Mai 1956.) 


Empfindlichkeit eisbildender Impfstoffe auf ihre Wasserlislichkeit 


Von RayMuUND SANGER, Ziirich!) 


Auf Grund von Messungen iiber die Eiskeimfahigkeit, ausgefiihrt an einer 
Mehrzahl von Stoffen mit Hilfe der Schaeferschen Wolkendiffusionskammer [1]?) 
hatten kiirzlich B. J. Mason und J. Hatier [2], [3] die Auffassung vertreten, 
dass entgegen den bis heute veréffentlichten Befunden innerhalb des Temperatur 
bereiches von 0 bis — 20°C in Wirklichkeit nur 2 Stoffe als eisbildungswirksam 
sich erweisen. Die beiden Stoffe umfassen die Mineralien Silberjodid (AgI) und 
Bleijodid (PbI,); alle die iibrigen bis heute als keimfahig gefundenen Stoffe 
verloren ihre Aktivitaét, wenn sie so weit gereinigt werden, dass auf spektrosko 
pischem Wege keine Verunreinigungen, insbesondere von Silber und Jod, nach 
gewiesen werden kénnen; dies trifft nach Mason unter anderem fiir: Kupfer- 
jodid (Cul) zu, das auch in unserem Laboratorium [4], [5] als keimfahig beobach 


1) Physikalisches Institut der ETH. 
*) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 540. 
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tet wurde. Ferner ziehen die beiden Verfasser nebenbei 
in Zweifel, dass das Mineral Kupfersulfid (CuS) die gleiche ausgepragte Eisbil- 
dungsfahigkeit wie AgI besitze, wobei sie allerdings darauf hinweisen, dass auch 
bei ihren Messungen CuS bei — 16° Anzeichen von Keimfahigkeit zeigte, doch 
k6nnte ihrer Meinung nach dies damit in Zusammenhang stehen, dass das ver- 
wendete CuS nur von der Reinheit der iiblichen kauflichen Stoffe war. Nach ihren 
Vorstellungen ist damit die bis anhin an andern Stoffen ausser AgI und PbI, 
uber — 20° beobachtete Gefrierkeimwirksamkeit, insbesondere diejenige von Cul 
und CuS, nur darauf zuriickzufiihren, dass diese Impfstoffe Verunreinigungen von 
Ag und I enthalten, und so iiber die Bildung des in kleinsten Spuren noch aktiven 
Agi Eiskeimfahigkeit vortauschen. Zugegebenermassen besteht bei Cul die Gefahr, 
dass auch bei sorgfaltiger Darstellung der Stoffe noch Spuren von AgI vorhanden 
sein kénnten. Nicht ohne weiteres ware es aber verstandlich, wie nach den 
ublichen chemischen Verfahren sorgfaltig dargestelltes CuS durch I verunreinigt 
sein kénnte. Tatsachlich haben nachtraglich vom analytisch-chemischen Labora- 
torium der ETH durchgefiihrte spurenanalytische Untersuchungen ergeben, dass in 
dem von uns verwendeten CuS keine Spuren von I (sicher kleiner als 1,0 - 10-® g I 
per Probe) nachzuweisen waren. 

Inzwischen hat dasselbe Laboratorium besondere Proben von CuS und Cul 
dargestellt; diese besitzen beziiglich ihres Metallgehaltes eine Reinheit, dass auch 
nach einem von O. GUBELI entwickelten kolorimetrischen Verfahren’) keine Ag- 
Spuren mehr aufzufinden und daher spektroskopisch AgI-frei anzusprechen sind. 
Eine erste behelfsmassige Priifung dieser Stoffproben auf ihre eisbildende Wirk- 
samkeit in einer gewohnlichen Wolkenkaltekammer haben gegensatzlich zu den 
Beobachtungen Mason und HALLETS ergeben, dass beide Impfstoffe gefrierkeim- 
aktiv sind, wobei die Einsatztemperatur mindestens fiir CuS in derselben Héhe wie 
bei Ag! liegt. Genauere Angaben iiber die Einsatztemperatur und tiber das Ausmass 
der Eisbildungswirksamkeit kénnen aber erst nach Fertigstellung der neuen Mess- 
anlage fiir die Ermittlung der Eiskeimfahigkeit von Iinpfstoffen gemacht werden. 

Unserer Auffassung nach besteht die Méglichkeit, dass die zur Schau getretene 
Zwiespialtigkeit der Messergebnisse der verschiedenen Beobachtergruppen physi- 
kalisch wohl begriindet ist. Es ist nicht gesagt, dass die Untersuchungen mit 
Hilfe der Schaeferschen Diffusionswolkenkammer zu denselben Ergebnissen 
fiihren wie jene mit Hilfe der gew6hnlichen Vonnegut-Schaeferschen Wolken- 
kAltekammer [6], namlich dann, wenn die untersuchten Impfstoffe eine zwar noch 
kleine, aber immerhin schon merkliche Wasserléslichkeit besitzen. 

In der Diffusionswolkenkammer 4ussern sich die von aussen eingeimptten 
Stoffpartikel unvermeidlich zuerst als Kondeusationskerne, bevor sie allfallig als 
Gefrierkerne zur Wirkung gelangen; dies auch dann, wenn die nee cheh, ae 
besonders ausgepragte Hygroskopizitat aufweisen. In den pelesrapece, und stets 
wachsenden Wassertrépfchen werden die eingebetteten Impfpartike , wenn ve 
Wasserléslichkeit zeigen, fortwahrend im Abbau begriffen sein und keine wo A 
definierte Oberflache mehr besitzen. Nach neueren Vorstellungen [4], U5] ae 
aber gerade der strukturelle atomare Aufbau der pbewerne ee agi 2 es 
eine iiberragende Rolle in der Frage der Eisbildungswirksamkeit; dieser ufbau 
wird aber durch die in den Tropfen sich abwickelnden Losungsvorgange eee 
lich gestért. Es ist daher nicht zu erwarten, dass Impfpartikel — pia ah Bs as 
vom Wasser unbenetzte Schwebeteilchen unterktihiten Wassertrop cn pee : ie 
den Gefriervorgang auslésen und daher im tiblichen Sinne gefrierkeimaktiv 


auch unsere Wahrnehmung 


Assi gent nek tralanalyse. 
3) Empfindlichkeit grdssenordnungsmassig 10fach gegeniiber der Spektralanaly 
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dasselbe auch tun, wenn sie in Wassertropfen eingelagert sind und ihre Ober-_ 
flache infolge ihrer Léslichkeit einem Abbauvorgang unterzogen wird. Solche 
Stoffe miissen daher, in die Diffusionswolkenkammer oder in die gewohnliche © 
Wolkenkaltekammer eingeimpft, nicht unbedingt zu den namlichen Ergebnissen ~ 
hinsichtlich ihrer Eisbildungswirksamkeit fiihren. Da AgI (Loslichkeit 10-7 g jem 
100 g H,O bei Zimmertemperatur) praktisch vollstandig wasserunléslich ist un 
Cul demgegeniiber (Léslichkeit 0,8 - 10~*/100 g H,O bei Zimmertemperatur) ver- 
haltnismAssig eine gute Léslichkeit zeigt, ware es damit verstandlich, wenn im erste= 
ren Falle beide Messverfahren dieselben Ergebnisse iiber die Gefrierkeimfahigkeit — 
zeitigen, wahrend im zweiten Fall widersprechende Befunde zutage treten. Auch © 
das Verhalten von CuS, das eine Léslichkeit von 3- 10~5 g/100 gH,O besitzt, a 
erklarlich. Denn das Wechselspiel der geschilderten, in der Oberflache der Impf-_ 
stoffe sich abwickelnden Vorgange muss sich notwendigerweise sehr mannnig-— 
faltig auswirken, so dass jede Art von Zwischenstufe sich einstellen kann, © 
wobei sicherlich auch das Ausmass der Hygroskopizitat mit ins Gewicht fallen — 
wird. » 
Das letzte Wort iiber die Eisbildungsfahigkeit ist noch nicht gesprochen; das _ 
Beobachtungsmaterial noch viel zu sparlich, um den Mechanismus der Wirk-_ 
samkeit eisbildender Kerne genau zu verstehen. So verfolgt diese Notiz nur 
den Zweck, darauf hinzuweisen, dass namhafte Griinde fiir die Zwiespaltigkeit der — 
Beobachtungen iiber die Eisbildungswirksamkeit der Impfstoffe in der Diffusions-_ 
wolkenkammer und in der gewohnlichen Wolkenkaltekammer geltend gemacht ~ 
werden kénnen und dass dabei die Wasserloslichkeit der untersuchten Impfstoffe — 
allfallig eine massgebende Rolle spielt. Wir stellen fest, dass bei den Unter-_ 
suchungen in den Wolkenkammern zu unterscheiden ist, ob sich die eisbildenden ~ 
Kerne den unterkiihlten Wassertropfen von aussen annahern oder ob sie sich 
bereits im Innern der Tropfen befinden. 
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Summary 


The ice-nucleating power of CuS and Cul, which H. PRuPpPACHER and the 
author discovered earlier, was confirmed with samples entirely free from Ag and I, 
and Ag respectively. It is stated that the same threshold temperatures are 
not necessarily found with a diffusion cloud chamber compared with a normal cloud 
chamber. The former displays the disadvantage against the latter since the nucleat- 
ing substance first acts as a condensation nucleus. The surface structure of the 


nuclei which determines the ice-nucleating power is more or less altered according 
to solubility. : 


(Eingegangen: 24. August 1956.) 
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} Die nachste Tagung der S.P.G. (Schweizerische Physikalische Gesellschaft) 
wird am 4,/5. Mai 1957 in Brunnen am Vierwaldstattersee stattfinden. Sekretariat: 
Basel, Klingelbergstrasse 82. K. P. MEYER 
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The Kinematics of Vorticity. Par C. TRuEsDELL (Indiana University Press, 
Bloomington 1954). XI + 232 pages, 5 figures; $6.00. 

Depuis quelques années, les travaux de TRUESDELL s’accumulent & un 
rythme imposant, et l’on peut bien dire que cet auteur a passé au rang des spécia- 
listes de l’hydrodynamique théorique. Son érudition est considérable, son style 
limpide et agréable. 

Le présent ouvrage rappelle en un certain sens la célébre Mécanique Analy- 
tique de WHITTAKER, en ce que l’ensemble des travaux originaux de quelque 
intérét sur la cinématique des fluides a été passé en revue et chaque auteur cité 
a sa place dans le développement de cette discipline. Cependant TRUESDELL ne 
veut pas faire un traité, mais un exposé original sur la formation, l’évolution et 
Yentretien des tourbillons. Nous ne saurions qu’approuver une telle entreprise, 
ayant nous-méme souligné son intérét dans un ouvrage sur les corps déformables 
paru il y a une dizaine d’années. TRUESDELL reprend les mémes themes princi- 
paux, tout en y mettant un soin et en y recherchant un détail dont il faut lui 
savoir gré; et comme nous l’avions esquissé dans ledit ouvrage, il développe tout 
Vappareil géométrique et analytique dont il a besoin, définit clairement la notion 
de «vorticité»!) et des mesures qu’il convient de lui attribuer, et développe tout 
au long la théorie de son évolution (en particulier: théorémes de BERNOULLI dans 
toute leur généralité), de sa convection et de sa diffusion. 

Ce n’est pas le cOté dynamique du mouvement dans les fluides qui retient son 
attention, mais le cété purement cinématique, aussi est-il clair que le tourbillon 
en est le centre d’attraction. 

La méthode de calcul employée est typique pour un ouvrage américain: pas 
un indice, pas une composante, bien qu’on nage dans les champs vectoriels et 
tensoriels. En revanche, emploi du calcul dyadique, que nous avons également 
développé et employé dans plusieurs ouvrages. Serait-ce peut-étre l’excellente 
Introduction to theoretical physics de LeicH Pacer, dans laquelle ont appris bien 
des physiciens théoriciens il y a 20 ou 30 ans, qui est responsable de cet usage ? 
A notre avis, c’est un usage excellent, qui permet de formuler les équations d’une 
facon trés condensée tout en évitant la fatigue que présentent les indices et les 
sommations. I] serait désirable que les auteurs européens se mettent a cette 


notation. 


1) Il vaudrait la peine de créer ce mot en frangais pour caractériser cette propriété des 
fluides de porter le champ du tourbillon. 
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Ce livre n’est pas fait pour des commengants. Il servira admirablement a | 
deux fins: aux spécialistes qui y trouveront une mine de renseignements bien | 
ordonnés, et aux professeurs ayant entre autres taches d’enseigner Vhydrodyna 
mique, car ils y trouveront des suggestions excellentes quant au choix des sujets” 
et 4 la forme de l’exposé. : 

L’ouvrage nous révéle l’importance des travaux ressortissant a de nombreuses — 
traditions: anglaise, francaise, italienne, allemande, russe, suisse, roumaine et | 
méme norvégienne. La typographie en est bonne et fait honneur aux Presses — 
universitaires d’Indiana qui, comme on le sait, impriment le Journal Rational — 
Mech. Anal., dont TRUESDELL est l’un des rédacteurs distingués. A. Mercier 


Materials for Nuclear Power Reactors. Von H.H.HAusNER und S.Bs 
RogorrF (Reinhold Publishing Corp., New York 1955). 224 S., 15 Fig.; $3.50. | 

Die fiir Leistungs-Kernreaktoren massgebenden physikalischen Grundlagen — 
sind heute so weit bekannt, dass sie nicht mehr das Hauptproblem bei der Weiter-_ 
entwicklung darstellen. Vielmehr sind es zur Zeit grundlegende Ingenieur- und 
Materialprobleme, welche in den Vordergrund geriickt sind und von deren Lésung 
der weitere Fortschritt abhangt. Das vorliegende Biichlein fiihrt auf sehr leicht-_ 
fassliche Art in die Werkstoffprobleme des Reaktorbaues ein. Ausgehend von 
einer kurzen Einleitung iiber die kernphysikalischen Grundlagen und die Kon-_ 
struktion von Kernreaktoren, werden die einzelnen Reaktorelemente und die fie 
sie in Frage kommenden Werkstoffe besprochen. Die Eigenschaften aller ver-_ 
wendungsfahigen Metalle und iibrigen Stoffe, wie zum Beispiel Graphit und 
Wasser, werden aufgeteilt nach ihrem Anwendungszweck (als Brennstoff, Mode-_ 
rator- und Reflektormaterial, Kiihlmittel, fiir Steuerung und Strahlenschutz) 
aufgefiihrt und an Hand eingehender Tabellen naher spezifiziert. Im ween 
sind die vor dem Genfer Kongress bekanntgewordenen Daten zusammengestellt. 
Eine Tabelle der chemischen Elemente und ihrer Isotopen, welche die Wirkungs- — 
querschnitte gegeniiber langsamen Neutronen angibt, sowie ein Glossary der 
im Reaktorbau tiblichen Fachausdriicke machen das Buch zu einem wertvollen 
Einfithrungswerk fiir alle jene, welche sich in materialtechnischer Hinsicht einen - 
Uberblick verschaffen méchten. W. Epprecht 


Servomechanisms and Regulating System Design, Bd. II. Von Ha- 
ROLD CHESTNUT und RoBEeRT W. Mayer (John Wiley & Sons, Inc., New York 
1955). 384 S., 217 Fig.; $8.50. 

Der vorliegende 2. Band erscheint in betrachtlichem Zeitabstand vom ersten. 
Der Grund lag darin, dass die Entwicklung, insbesondere der nichtlinearen Servo-_ 
technik, sehr rasche Fortschritte erfahren hat, so dass ein Teil der bereits im 
Manuskript vorliegenden Kapitel vollstandig neu bearbeitet werden mussten. 
Der 2. Band wendet sich an den Konstrukteur von Servosystemen, wobei viele 
Einzelprobleme eingehender behandelt werden. Ein besonderer Abschnitt, nam 
lich das 1. Kapitel, widmet sich der Messtechnik von Teilen oder auch ganzen 
Servomechanismen, wobei sowohl die zeitlichen Einschwingvorgange wie auch 
der Frequenzgang untersucht werden. Im 2. Kapitel wird der Einfluss der For m 
des Eingangssignals auf den Entwurf eines Servosystems diskutiert. Hierzu_ 
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Geschwindigkeit und Leistung erfiillt werden kénnen. Im 4. Kapitel wird eine 
Ubersicht von Korrekturnetzwerken gegeben, welche eine bestimmte Frequenz- 
abhangigkeit der Verstarkung aufweisen. Im 5, Kapitel wird der Entwurf des 
Verstarkers diskutiert. Besondere Aufmerksamkeit wird der Nullpunktkonstanz 
von Gleichstromverstarkern gewidmet. In diesem Zusammenhang sind auch die 
Tragersysteme wichtig, welche die Gleichstromgrésse in eine Wechselstromgrisse 
umsetzen, so dass diese verstirkt und wieder in eine Gleichstromgrésse zuriick- 
verwandelt werden kann. Eine besondere Rolle spielen die reinen Wechselstrom- 
servomechanismen, welchen das 6. Kapitel gewidmet ist. Die Stabilisierung und 
der Einfluss von Abweichungen der Tragerfrequenz werden eingehend behandelt. 

Wahrend bis hierher rein lineare Vorgange betrachtet wurden, sind die fol- 
genden Kapitel nichtlinearen Vorgangen gewidmet. So wird im 7. Kapitel die 
Linearisierung von nichtlinearen Elementen fiir kleine Abweichungen diskutiert. 
Diese Methode wird fiir grosse Ausschlage nicht mehr Giiltigkeit haben, so dass 
fiir solche Vorgange eine andere Linearisierung notwendig sein wird, was seinen 
Niederschlag im 8. Kapitel findet. Wichtig sind hier die Sattigungserscheinungen 
und die Diskussion der Stabilitatsfragen. Die Sattigungseffekte kénnen bei Posi- 
tions-Servosystemen einen wesentlich andern Einschwingvorgang zur Folge haben, 
als dies bei linearen Systemen der Fall ware. Besonders wichtig ist in diesem 
Zusammenhang auch das Spiel, das sich durch Reibung und Ungenauigkeiten der 
Getriebe ergibt. Bereits in diese Kategorie gehéren die Relais-Servosysteme, die 
einen besonderen Fall der nichtlinearen Systeme darstellen. Im letzten Kapitel 
wird die Anwendung von nichtlinearen Elementen in Servosystemen behandelt, 
und es wird gezeigt, dass mit solchen Elementen unter Umstanden eine bedeu- 
tende Verbesserung des Verhaltens erzielt werden kann. Eine Bibliographie bis 
1954, auf welche im Text Bezug genommen wird, erganzt das Buch. 

Der Vorzug des Buches besteht darin, dass eine ausgezeichnete Auswahl der 
Probleme getroffen wurde, um die wesentlichsten Gesichtspunkte herauszuar- 
beiten. Das Buch kann jedem, der sich eingehender mit dem Entwurf von Servo- 
‘systemen abgeben muss und der die grundlegenden Aspekte des ersten Bandes 
kennt, als Wegweiser wertvolle Dienste leisten. H. Weber 


Reaktionen in und an festen Stoffen. Von K. Haurre (Springer-Verlag, 
Berlin 1955). 696 S., 427 Abb.; DM 78.-. 

Die Halbleitertheorie hat sich in den letzten Jahren als neues Bindeglied 
zwischen Physik und anorganische Chemie eingefiigt. Die chemischen Reaktionen 
im festen Zustand kénnen heute auf Grund der Fehlordnungs- und Halbleiter- 
phanomene sehr weitgehend verstanden und in sehr vielen Fallen bereits auch 
quantitativ erfasst werden. Es ist daher sehr zu begriissen, dass K. HAvFFE den 
Versuch gewagt hat, die Fehlordnungs- und Halbleitertheorien als Grundlage fiir 
die Erérterung der Reaktionen an und in festen Korpern zu wahlen. Sein Werk 
erfiillt dadurch einen doppelten Zweck: es zeigt dem Chemiker, wie das chemische 
Geschehen im festen Zustand studiert werden kann und erinnert anderseits den 
Physiker daran, dass die Halbleitertheorie nicht nur fiir die heute in den Vorder- 
grund geriickten Transistoren und Gleichrichter nutzbringend ist, sondern eine 
viel weiter gespannte Bedeutung besitzt. Sie reaktivierte die eine Zeitlang in den 

‘Hintergrund geriickte Festkérperchemie sehr stark. Das Buch dart schon heute 
als Standardwerk bezeichnet werden, das jedem Festkorperchemiker und -phy- 
siker wertvolle Dienste zu leisten vermag. 

Der erste Teil fiihrt zunaichst in die Fehlordnungserscheinungen von Tonen- 

‘und Valenzkristallen ein und gibt die Berechnungsgrundlagen wie auch die von 


» 


Press 1955). 265 pp., 61 figs.;°35 s. ; 


_ which he is unfamiliar. Indeed, in the short time that this reviewer has had 
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ScHotrKy und WaGNER eingefiihrte Symbolik an. Sodann werden an Hand gut. 
untersuchter Beispiele die Beziehungen zwischen elektrischer Leitfahigkeit sowie 
der Thermokraft und den Fehlstellen dargelegt und die Versuchsmethoden, 
welche zur Erforschung dieser Beziehungen tiblich sind, erlautert. Als erste An-— 
wendung wird das Problem der Randschichterscheinungen besprochen, das beson-— 
ders fiir die Deutung der Chemisorption und damit fiir die Katalyseforschung von” 
Interesse ist. Ein sehr eingehendes Kapitel setzt sich anschliessend mit der 
Diffusion auseinander, wobei wie in den vorangehenden Abschnitten sehr viele 
Beispiele zur Darlegung herangezogen werden. Im zweiten Teil des Werkes_ 
werden die gewonnenen Erkenntnisse auf die in der Praxis bedeutungsvollen 
chemischen Reaktionen von Festkérpern angewandt, so die Oxydationsvorgange © 
-an metallischen Stoffen, die Bildung von Ionenverbindungen hoéherer Ordnung 
(zum Beispiel Spinellbildung) und schliesslich fiir die Reduktions- und Rést- 
prozesse. Dadurch, dass hier wie schon in den grundlegenden Kapiteln ausser- 
ordentlich viele Beispiele und Literaturangaben verarbeitet und durch gute ~ 
Figuren veranschaulicht sind, wird das Buch zu einem wertvollen Nachschlage-— 
werk, was besonders auch durch die gut redigierten Autoren- und Sachregister © 
betont wird. Das Werk bleibt indessen ein sehr gutes Lehrbuch, enthalten doch ~ 
alle Kapitel einfiithrende und zugleich zusammenfassende Abschnitte, die zusam-_ 
mengenommen einen vorziiglichen Uberblick iiber die in den speziellen Abschnit- — 
ten durch Beispiele erhdrteten Forschungsergebnisse geben. Trotz der Fiille des — 
Stoffes bleibt das Werk iibersichtlich und leichtfasslich. So schliesst es denn 
eine bisherige Liicke aufs beste, und es ist daher jedem zu empfehlen, der sich 
- mit Halbleiterproblemen, Reaktionen im festen Zustand oder Wechselwirkungen ~ 
zwischen Festkoérpern und Gasen zu befassen hat, sich in dieses Werk zu vertiefen. — 
W. Epprecht 


Nuclear Magnetic Resonance. By E. R. ANDREw (Cambridge University 


Professor ANDREW has written an excellent review of the subject of nuclear 
magnetic resonance. It is particularly suited for those who wish to begin their 

. study of this field, since the book is not only comprehensive, but the physical - 
ideas have been explained with great clarity and with a minimum of mathe- — 
matics. For details, the reader is referred to the original papers; Professor 
ANDREW’s listing of these papers forms a most comprehensive bibliography 
which in itself is a valuable contribution to the book. 7 
The book, however, is not only suited for beginners, but for practicing 
physicists as well, since one often desires a quick review ina part of the field in — 


the book in his possession, he has found Professor ANDREW’s book exceptionally — 
_ well suited to that purpose, since, in addition to the bibliography, a most compre- _ 
hensive listing of nuclei and sample materials (organic and inorganic solids, 
metals) which have been investigated by nuclear magnetic resonance is included 
in the appendices. Bee 
The only drawback lay outside the author’s hands: it is always difficult to 
write a book in a rapidly developing field. As a consequence, several of the 
newest and most exciting developments in the field are not covered. Professor 
ANDREw, however, is aware of this inevitable shortcoming, bringing, for example, 
by means of a footnote written as the book went to press, the reader’s attention _ 


_ to the results of studies at extremely high resolution, typified by the work 
ARNOLD and ANDERSON. W. G. Proctor 
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Der Schatzkanzler auf der Darius-Vase zahlt die eingenommenen Steuern auf dem Rechenbrett 
zusammen. Aus B. L. VAN DER WAERDEN, Erwachende Wissenschaft 


Erwachende Wissenschaft 


Agyptische, babylonische und griechische Mathematik © 


Von Dr. B. L. VAN DER WAERDEN, Professor an der Universitat Zurich. Aus dem 

Hollandischen iibersetzt von Hetca Hasicut und mit Zusdtzen vom Autor ver- 

sehen, 488 Seiten mit 180 Figuren. In Ganzleinen gebunden Fr./DM 37.50. Samm- 
lung « Wissenschaft und Kultur», Band 8. 

Aus dem Inhalt: Einfiihrung: Wozu Geschichte der Mathematik — 1. Die Agypter - 

Il. Zahlensysteme, Ziffern und Rechenkunst — III. Babylonische Mathematik — 

IV. Das Jahrhundert von Thales und Pythagoras — V. Das goldene Zeitalier — V1. Die 


. Zeit Platons — VII. Die Alexandrinische Zeit — VIII. Niedergang der griechischen 


Mathematik — Sachverzeichnis, 


Der Autor war bestrebt, das Buch einerseits wissenschaftlich exakt, andererseits 
aber auch leicht fasslich zu machen in dem Sinne, dass jeder, der an einer Mittel- 
schule Mathematik gelernt hat und sich fiir die Geschichte der Mathematik inter- 
essiert, es verstehen kann. — Die Geschichte der Mathematik erscheint im rech- 
ten Licht nur, wenn man sie gegen den Hintergrund der allgemeinen Kultur- 


geschichte betrachtet. Das Leben in den griechischen Stadten in der klassischen — 


Zeit ist anders als das Leben an den Fiirstenhéfen von Alexandrien und Syrakus, 
und die Wissenschaft wird vom Leben manchmal entscheidend beeinflusst. 
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